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ABSTRACT 

Magnetization  measurements  were  carried  out  in  a study  of 
local  magnetic  effects  (magnetic  clusters,  paramagnons,  etc.) 
in  the  transition-group  alloy  systems,  Ni-Rh,  Pd-Ni,  and  Fe-Cr, 
near  their  critical  compositions  or  critical  temperatures  for 
ferromagnetic  ordering.  Similar  measurements  were  made  of  the 
detailed  magnetization  processes  associated  with  the  micto- 
magnetic  (spin-glass)  ordering  of  Cu-Mn  alloys.  Also  conducted 
was  a magnetic  and  neutron  diffraction  investigation  of  high- 
order  (biquadratic)  exchange  interactions  in  the  rare-earth 
compounds,  PrAg  and  DySb. 
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INTRODUCTION 
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The  magnetism  of  alloys  containing  metals  of  the  first  trans- 
ition group  (the  iron  group)  is  a major  unresolved  area  of  solid 
state  physics  in  which  there  is  considerable  current  interest. 
Although  definitive  advances  have  been  made  towards  an  under- 
standing of  isolated  magnetic  impurity  atoms  in  a metallic  host, 
the  more  complex  problem  concerning  the  magnetic  states  of  con- 
centrated, atomically  disordered  alloys  continues  to  pose  important 
basic  questions.  Some  of  these  questions,  however,  appear  to  be 
answered  by  the  recent  discoveries  of  intrinsically  inhomogeneous 
magnetic  states  in  certain  concentrated  nickel-based  alloys  (spe- 
cifically, Ni-Cu'*'  ^ Ni-V,^'"^  and  Ni-Rh^  . A new  phenomenological 
picture  has  emerged  for  the  magnetism  of  these  alloys,  in  which  the 

formation  and  size  of  a stable  magnetic  moment  of  a transition- 

11-14 

group  atom  are  critically  determined  by  the  local  environment. 

The  broad  statistical  distribution  of  local  environments  in  these 
solid-solution  alloys  gives  rise  to  ferromagnetic  and  paramagnetic 
states  that  are  spatially  inhomogeneous.  However,  much  remains 
to  be  learned  about  the  various  local  conditions  for  stable  or 
nearly  stable  magnetic  moments  in  different  systems. 

In  a related  class  of  transition-group  alloys  known  as  spin 
glasses  or  mictomagnets , the  ordered  magnetic  state  is  also 
thought  to  be  spatially  inhomogeneous.  However,  the  inhomogeneity 
in  this  case  refers  mainly  to  the  directions  rather  than  the 
magnitudes  of  the  moments  of  the  randomly  situated  magnetic  atoms. 
This  orientational  inhomogeneity  is  considered  to  arise  from  an 


A 


oscillation  in  the  sign  of  the  exchange  coupling,  between 
ferromagnetic  and  antiferromagnetic,  as  a function  of  distance. 

For  the  prototypal  spin-glass  systems,  Cu-Mn  and  Au-Fe,  it  has 
recently  been  observed^ ' ^ that  the  very  pronounced  irreversi- 
bilities in  their  magnetic  properties  (known  from  earlier 
18  19 

work)  ' set  in  precisely  at  the  so-called  freezing  temperature, 

20 

which  is  marked  by  a sharp  peak  in  the  initial  susceptibility 

and  depends  on  the  alloy  composition.  The  irreversibility  of  the 

magnetization  with  field  and  temperature  has  been  tentatively 

18  21  2~> 

attributed  to  some  type  of  local  magnetic  anisotropy,  ' ' 

but  the  underlying  microscopic  mechanism  is  still  unknown.  At 
present,  there  are  conflicting  opinions  even  about  the  nature  of 
the  mictomagnetic  ( ferro-antiferromagnetic)  ordering  or  freezing 
process.  Thus,  many  of  the  basic  questions  about  the  magnetism 
of  these  alloys  remain  unresolved. 

With  regard  to  transition-group  alloys,  our  research  program 
involved  detailed  experimental  studies  of  the  thermodynamic  prop- 
erties of  Ni-Rh,  Pd-Ni,  Fe-Cr,  and  Cu-Mn  alloys  near  their 
critical  compositions  or  temperatures  for  the  onset  of  ferro- 
magnetic or  mictomagnetic  order.  The  summary  purpose  of  these 
investigations  was  the  determination  of  the  essential  inhomo- 
geneous character  of  the  ordered  magnetic  states  and  the 
paramagnetic  states  from  which  they  derive  in  concentrated  solid- 
solution  alloys. 

The  magnetisi..  : rare-earth  metals  and  intermetallic  compounds, 
compared  to  that  of  their  transition-group  counterparts,  is  quite 
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well  understood.  Almost  without  exception,  the  magnetic  moments 

of  rare-earth  atoms  in  a metal  derive  from  well-localized  4f 

electron  states  and  are  coupled  together  via  indirect  exchange 

interactions  involving  the  itinerant  conduction  electrons. 

Furthermore,  even  in  a metal,  the  crystal-field  energy  levels 

for  a rare-earth  atom  can  readily  be  determined  from  inelastic 

neutron  scattering  experiments  and,  consequently,  the  magnetic 

polarization  of  the  atom  can  be  calculated  as  a function  of  the 

total  effective  (net  exchange  plus  external)  field,  the  temperature, 

and  its  orientation  in  the  crystal.  Such  a procedure  was  followed 

in  a recent  study  of  a polycrystalline  sample  of  PrAg,  and  the 

paramagnetic  data  revealed  that  the  net  exchange  coupling  in  this 

intermetallic  compound  contains  a higher-order  (biquadratic)  com- 

23 

ponent  comparable  in  strength  to  the  normal  bilinear  component. 

Our  magnetic  and  neutron  diffraction  experiments  have  shown 
that  the  anomalous  magnetic  properties  of  PrAg  below  as  well  as 


above  its  Neel  point  can  be  directly  attributed  to  the  existence 


of  biquadratic  interactions. 
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Our  magnetic  studies  have  also 


revealed  that  strong  biquadratic  coupling  exists  in  the  compound, 
DySb.25'26 
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EXPERIMENTAL  RESULTS  AND  DISCUSSION 


Transition-Group  Alloys 

During  the  contract  period,  our  work  pertinent  to  transition- 
group  alloys  consisted  of  the  following:  a detailed  magnetization 
study  of  the  Ni-Rh  system  near  the  critical  composition  for  ferro- 
magnetism and  a similar  but  preliminary  study  of  the  Ni-Pd  system; 
a study  of  the  magnetic  properties  of  Fe-Cr  alloys  near  their 
critical  temperatures;  magnetization  measurements  on  a Cu-Mn 
alloy  near  and  below  the  freezing  temperature  for  its  micto- 
magnetic  state. 

Our  magnetization  studies  of  Ni-Rh  began  about  five  years 

ago,  and  the  preliminary  results  of  this  work  have  already  been 
8 9 

reported.  ' This  investigation  has  subsequently  been  pursued 
in  greater  detail  over  a much  broader  range  of  alloy  compositions, 
extending  through  the  critical  composition  for  ferromagnetism. 

The  data  analysis  for  this  detailed  investigation  was  completed 
about  a year  ago,  and  a full  report  of  this  work  was  recently 
published  in  the  Physical  Review;^®  a reprint  copy  of  this  paper 
is  attached  as  Appendix  A. 

In  summary,  our  Ni-Rh  results  show  that  the  paramagnetic  sus- 
ceptibility can  be  resolved  into  a Curie-Weiss  component  and  a 
weakly  temperature  dependent  component  of  the  exchange-enhanced 
Pauli  type.  The  onset  of  ferromagnetism  is  found  to  be  directly 
related  to  a divergence  of  the  Curie-Weiss  susceptibility  com- 
ponent and  is  attributed  to  interactions  between  superparamagnetic 
clusters.  The  average  magnetic  moment  per  cluster  is  deduced  to 
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be  ^20-24  over  the  composition  range  investigated  (51-65  at.%Ni). 
The  concentration  of  the  magnetic  clusters,  though  always  dilute, 
increases  rapidly  as  the  critical  composition  is  approached  from 
the  paramagnetic  side.  Moreover,  the  cluster  concentrations  cor- 
relate closely  with  the  statistical  occurrences  of  extremely 
Ni-rich  local  regions  in  these  atomically  disordered  alloys.  The 
intrinsically  inhomogeneous  paramagnetic  and  ferromagnetic  states 

of  Ni-Rh,  thus  described,  are  basically  similar  to  those  deduced 

4 6 7 

earlier  for  Ni-Cu  and  Ni-V  ' alloys  near  their  critical  compo- 
sitions. These  results  disagree  with  previous  speculation  that 

the  exchange  enhancement  in  Ni-Rh  does  not  become  critical  until 

27 

ferromagnetic  ordering  takes  place.  However,  our  evidence  for 
the  existence  of  giant-moment  clusters  in  paramagnetic  Ni-Rh, 
which  presumably  arise  from  critical  local  enhancement,  is  con- 
sistent with  recent  heat  capacity  measurements  in  a magnetic 
field. 28 

Paramagnetic  Ni-Pd  alloys  have  been  thought  to  provide  an 
ideal  isoelectronic  situation  for  the  occurrence  of  paramagnons 

(spin  fluctuations)  associated  with  a strong  exchange  enhancement 
29 

at  the  Ni  sites.  Indeed,  the  results  of  our  preliminary 
magnetization  study  concur  with  this  expectation  for  Ni-Pd  alloys 
very  dilute  in  Ni.  However,  as  the  critical  composition  for 
ferromagnetism  is  approached  more  closely,  our  results  reveal  the 
appearance  of  superparamagnetic  clusters,  which  signifies  that 
the  local  exchange  enhancement  at  certain  Ni  atoms  attains  criti- 
cality even  though  the  alloy  remains  magnetically  disordered.  Our 


L. 
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preliminary  results  on  Ni-Pd  were  reported  at  the  1976  March 
Meeting  of  the  American  Physical  Society;  the  abstract  is  at- 
tached as  Appendix  B. 

We  have  measured  the  detailed  magnetization-field-temperature 
properties  of  several  Cr-rich  Fe-Cr  alloys  close  to  their  ferro- 
magnetic Curie  points.  Our  results  show  that  the  critical 

exponents  have  anomalous  values,  similar  to  those  reported  earlier 

30  9 

for  weakly  ferromagnetic  Pd-Fe  and  Ni-Rh  alloys.  A further 
and  more  revealing  similarity  between  the  latter  alloys  and  Fe-Cr 
was  found  in  the  values  of  the  critical  coefficients,  which  show 
that  weakly  ferromagnetic  Fe-Cr  alloys  also  have  their  magnetiza- 
tion distributed  inhomogeneous ly  in  giant  magnetic  clusters. 

The  results  of  this  work  were  presented  at  the  1976  International 

Conference  on  Magnetism;  a copy  of  the  paper  that  appeared  in  the 

31 

conference  proceedings  is  attached  as  Appendix  C. 

In  an  investigation  of  Cu-Mn,  we  have  explored  the  detailed 
magnetization-vs-f ield  behavior  of  a v20  at.%Mn  alloy  cooled  in 
zero  field  to  various  temperatures.  The  temperature  range  of 
measurement  extended  through  T , the  so-called  freezing  temperature 
for  the  mictomagnetic  state,  at  which  the  initial  susceptibility 
has  a sharp  maximum.  Our  results  show  Tm  is  also  marked  by  a 
qualitative  change  in  the  nature  of  the  magnetization  processes. 

Below  T , the  changes  in  magnetization  above  a certain  threshold 
field  become  irreversible,  exhibiting  hysteresis  and  remanence. 

The  value  of  this  threshold  field  and  the  amount  of  irreversibility 
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(i.e.,  the  area  of  the  hysteresis  loops)  both  decrease  gradually 
with  increasing  temperature  and  go  to  zero  at  T , above  which 
the  magnetization  is  reversible  at  all  fields.  In  fact,  it 
appears  that  the  susceptibility  peak  at  Tm  is  itself  a manifes- 
tation of  the  onset  of  magnetic  irreversibility,  which  is  also 

16 

suggested  by  recent  work  on  Au-Fe.  This  effect  is  quite  central 

and  must  be  taken  into  account  in  any  theoretical  model  for  a 

32 

spin-glass  or  mictomagnetic  alloy.  A brief  report  of  our  work 

17 

was  recently  published  in  Solid  State  Cimmunications ; a copy 
of  this  paper  is  attached  as  Appendix  D. 

Rare-Earth  Compounds 

During  the  contract  period,  we  have  shown  from  theoretical 
analysis  that  the  anomalous  magnetic  properties  of  the  cubic 
(CsCl-type)  compound  PrAg  can  be  explained  on  the  basis  of  a 
strong  biquadratic  exchange  coupling.  Furthermore,  we  have 
completed  a detailed  magnetization  study  of  the  NaCl-structured 
compound  DySb,  whose  magnetic  properties  suggest  a similar  pre- 
dominance of  higher-order  (presumably  biquadratic)  exchange. 

Our  active  interest  in  PrAg  began  shortly  after  neutron  dif- 
33 

fraction  experiments  had  shown  this  compound  to  be  an 
antiferromagnet , contrary  to  previous  reports  of  ferromagnetism, 
and  had  resulted  in  a determination  of  its  crystal-f iedl  energy 
levels.  We  embarked  on  a magnetization  study  of  a poly-crystal- 
line sample  of  PrAg  (since  single  crystals  were  then  unavailable) 


k. 
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and  our  results,  some  of  which  have  been  briefly  reported, 

can  be  summarized  as  follows:  PrAg  is  indeed  an  anti-ferro- 

magnet  with  a Neel  point  (T„)  of  vH°k.  At  all  temperatures 

below  TN,  it  undergoes  a spin-flop  transition  at  a critical 

field  (H  ) of  ^5  kOe.  Above  H , its  magnetization  (M)  rises 
c c 

approximately  as  the  cube  root  of  the  field  (H) , in  contrast 
with  the  normal  linear  behavior.  Moreover,  the  initial  suscep- 
tibility increases  anomalously  below  its  peak  at  T^.  For  the 
paramagnetic  state  above  TN,  our  M vs . H data  were  compared 

against  M vs.  (total  effective  field)  calculated  for  the 

33 

crystal- field  states  deduced  from  the  neutron  measurements. 
Within  the  effective  field  approximation,  He^  - H at  a given 
M represents  the  net  exchange  field  on  each  Pr  atom,  and 

would  normally  be  expected  to  vary  linearly  with  M.  Instead, 
we  found  that 


exch 


(1) 


with  A'  negative  and  of  the  same  order  of  magnitude  as  X. 

In  our  analysis  of  this  behavior,  it  was  recognized  that 

the  X'M^  term  in  Eq.  (1)  can  arise  from  a term  in  the  exchange 

Hamiltonian  that  is  quartic  in  the  spin  variables.  The  simplest 

-+  2 

such  term  has  the  isotropic  biquadratic  form  (S^  • S^)  . It  was 

therefore  assumed  that  the  exchange  Hamiltonian  of  the  ith  Pr 
atom  interacting  with  j Pr  neighbors  may  be  written  as 


Hi  " ‘ \ Jij  Si 


S3  - ^ 


(3, 


V2 


(2) 
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where  the  first  term  on  the  right  represents  the  usual  bilinear 

I 

exchange.  Applying  this  Hamiltonian  with  negative  J\  ^ to  the 
triplet  ground  state  of  PrAg  (the  excited  states  being  ignorable 
except  at  very  high  temperatures)  within  the  molecular  field 
approximation,  we  found  that  Eq.  (1)  is  totally  justified  for 
the  paramagnetic  state.  Furthermore,  this  model  Hamiltonian 
was  found  to  generate  all  the  anomalous  antiferromagnetic 
properties  of  PrAg  that  are  described  above.  A full  report  cn 
this  work,  prepared  in  collaboration  with  T.  0.  Brun  and  G.  H. 

Lander  of  the  Argonne  National  Laboratory,  was  recently  published 

24 

in  the  Physical  Review;  a copy  of  this  paper  is  attached  as 
Appendix  E. 

The  compound  DySb,  which  becomes  antiferromagnetic  at 
Tn  ^ 9°K  via  a first-order  transition,  34,35  undergoes  rapid  field- 
induced  transitions  whose  characteristics  depend  on  the  crystal- 
lographic direction  of  the  applied  field.  With  increasing 
field  along  <100>,  there  is  first  a transition  to  a state  whose 
magnetization  corresponds  to  v 5uB/Dy-atom  and  then  another  to  a 
state  of  v lOjig/Dy-atom,  corresponding  to  complete  ferromagnetic 
alignment.  The  magnetization  of  the  intermediate  state  suggests 
that  half  of  the  Dy"3"1”  moments  are  parallel  to  the  field  and  the 
other  half  are  orthogonal. 

From  magnetization  and  neutron  diffraction  studies  of  a 
DySb  crystal  (the  latter  study  conducted  at  the  Argonne  National 
Laboratory  in  collaboration  with  T.  0.  Brun  and  G.  H.  Lander), 
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we  have  determined  the  various  magnetic  states  of  this  material 
and  its  magnetic  phase  diagrams  for  fields  along  different  crys- 
tallographic directions.  In  particular,  we  have  found  that  in 
the  intermediate-field  ordered  state  of  DySb  the  magnetic  moments 
do  indeed  form  an  orthogonal  configuration,  which  can  be  attrib- 
uted to  strong  biquadratic  coupling  of  appropriate  sign,  i.e., 

I 

J — in  Eq.  (2)  is  negative,  as  in  PrAg.  The  results  of  this  work 

were  presented  at  the  1974  Conference  on  Magnetism  and  Magnetic 

25 

Materials  and  were  published  in  the  conference  proceedings;  a 

copy  of  this  paper  is  attached  as  Appendix  F. 

To  further  probe  the  existence  of  higher-order  interactions 

in  DySb,  we  pursued  our  detailed  magnetization  measurements 

throughout  the  paramagnetic  regions  of  field  and  temperature. 

The  paramagnetic  data  obtained,  when  analyzed  as  described  above 

for  PrAg,  showed  that  the  coefficient  of  the  effective  biquadratic 

coupling  in  DySb  depends  very  sensitively  on  the  temperature  and 

on  the  crystallographic  direction  of  the  field.  These  results, 

which  we  are  presently  studying  theoretically,  were  reported  at 

the  1976  International  Conference  on  Magnetism.  A copy  of  the 

2 6 

paper  that  appeared  in  the  conference  proceedings  is  attached 
as  Appendix  G. 
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Magnetic  properties  of  Ni-Rh  alloys  near  the  critical  composition  for  ferromagnetism* 

W C.  Muellner'  and  J.  S.  Kouvel 
Department  of  Physics.  University  of  Illinois.  Chicago.  Illinois  606 SO 
(Received  22  January  1975) 

The  magnetizations  of  Ni-Rh  alloys  on  either  side  of  the  critical  composition  for  ferromagnetism 
(c,.nt  % 63-at.%  Nil  were  measured  between  4.2  and  250  “k  in  fields  up  to  56  kOe.  The  initial 
paramagnetic  susceptibility  of  each  alloy  is  shown  to  be  resolvable  into  a Curie-Weiss  component  and  a 
weakly  temperature-dependent  component  of  the  exchange -enhanced  Pauli  type.  The  latter  component, 
which  is  essentially  equal  to  the  high-field  differential  susceptibility  at  low  temperatures,  is  maximum 
near  ccnl.  The  onset  of  ferromagnetism,  however,  is  more  directly  related  to  divergence  of  the 
Cune-Weiss  susceptibility  component  and  is  attributed  to  interactions  between  superparamagnettc 
clusters.  The  average  moment  per  cluster  is  deduced  to  be  — (20-24)/z„  over  the  composition  range 
investigated  (65-51-at.%  Ni)  The  concentration  of  the  magnetic  clusters,  though  remaining  dilute, 
increases  rapidly  as  c£tlI  is  approached  from  the  paramagnetic  side,  moreover,  it  correlates  closely 
with  the  statistical  occurrence  of  extremely  Nt-nch  local  regions  in  these  atomically  disordered  alloys. 

Comparisons  are  made  with  the  Nt-Cu  and  Ni-V  systems  near  c„„ . for  which  recent  magnetization 
studies  have  revealed  the  analogous  existence  of  dilute  concentrations  ol  giant  magnetic  clusters. 


INTRODUCTION 

The  alloys  of  nickel  pose  a special  challenge  to 
our  understanding  of  magnetism  in  metallic  sys- 
tems. Historically,  their  spontaneous  magnetiza- 
tions and  related  properties  were  once  thought  to 
have  found  a simple  explanation  in  the  collective- 
electron  rigid-band  model,  particularly  in  the  case 
of  nickel  alloys  with  non- transition-group  metals 
(Cu,  Zn,  etc. ).  1 This  theoretical  rationale  was 
applied  fairly  successfully  to  ferromagnetic  alloys 
all  the  way  out  to  the  critical  composition  (c£rit, 
where  the  Curie  point  Tc-0  °K)  and  even  beyond 
into  the  paramagnetic  composition  range,  where 
the  exchange  forces  are  expected  to  cause  an  en- 
hancement of  the  Pauli  susceptibility.  However, 
it  has  been  known  that  certain  detailed  properties 
of  the  prototypal  Ni-Cu  system  near  ccrl,  (e.g. , 
pronounced  susceptibility  increases  at  lew  tem- 
peratures,2,3  an  essentially  temperature-inde- 
pendent component  of  the  specific  heat*'5)  are  dis- 
tinctly anomalous  from  a rigid-band  model  view- 
point and  are  more  readily  attributable  to  super- 
paramagnetic  clusters  acting  as  giant  local  mo- 
ments.5 Generally,  these  anomalies  have  been 
dismissed  as  extraneous  manifestations  of  gross 
chemical  inhomogeneities  in  the  alloy  samples.7 
The  fact  that  Ni,  as  a dilute  impurity  in  Cu  or  any 
other  nonmagnetic  metal  host,  shows  no  signs  of 
having  a stable  local  moment  (unlike  Fe  or  Mn)8 
seemed  to  argue  against  the  possibility  of  local  mo- 
ments as  intrinsic  phenomena  in  concentrated  Ni 
alloys. 

Many  recent  developments  have  changed  this  pic- 
ture profoundly,  especially  with  regard  to  Ni-Cu. 
Most  recently,  neutron  scattering  experiments  on 
ferromagnetic  Ni-Cu  alloys  by  Aldred  ct  nl.*  (and 


earlier  on  one  alloy  by  Cable  ct  nl . :o)  have  demon- 
strated that  the  3z/-electron  moment  of  a Ni  atom 
depends  sensitively  on  its  local  environment,  de- 
creasing in  magnitude  with  increasing  number  of 
neighboring  Cu  atoms.  Indeed,  these  results3  in- 
dicate a gradual  evolution  with  increasing  Cu  con- 
centration toward  the  highly  inhomogeneous  spatial 
distributions  of  magnetization  in  the  weakly  ferro- 
magnetic alloys  near  Co.,,.  These  distributions,  as 
deduced  previously  by  Hicks  ct  nl.  from  similar 
experiments,  u consist  of  discrete  regions  of  mag- 
netization, called  “magnetic  polarization  cloucis,  ” 
each  extending  over  many  atoms  and  having  an 
average  total  moment  of  about  10ua.  Later  sus- 
ceptibility measurements  on  Ni-Cu  revealed  that 
the  giant  polarization  clouds  persist,  as  superpara- 
magnetic  entities,  well  into  the  paramagnetic  com- 
position range.  12  Moreover,  the  very  dilute  con- 
centrations of  these  superparamagnetic  clouds  in- 
dicated that  they  nucleate  at  local  regions  that  are 
extremely  Ni  rich.  These  local  Ni-rich  regions, 
the  cores  of  the  polarization  clouds,  may  occur 
statistically  even  in  perfectly  random  solid  solu- 
tions, as  has  been  discussed  by  Perrier  ct  nl.  13; 
the  short-range  atomic  clustering  that  normally 
exists  in  Ni-Cu  alloys9’1'  simply  increases  the  con- 
centration of  these  regions.  This  statistical  pic- 
ture, with  the  magnitude  of  each  Ni  atomic  moment 
having  a simple  prescribed  dependence  on  its  local 
chemical  environment,  gives  a good  quantitative 
fit  to  the  saturation  magnetizations  of  Ni-Cu  alloys 
near  cerlt.  13  An  analogous  but  more  complicated 
prescription  adopted  by  Robbins  ct  nl . 15  gives  ex- 
cellent agreement  over  the  entire  Ni-Cu  composi- 
tion range.  This  viewpoint  of  a spatially  inhomo- 
geneous magnetic  state  has  been  incorporated,  of- 
ten with  some  allowance  for  the  magnetic  aspects 
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ol  the  local  environments,  into  several  theoretical  EXPERIMENTAL  TECHNIQUES 

analyses  of  the  magnetization  and  related  proper- 
ties of  Ni-Cu.  10-23  Nickel-rhodium  alloy  buttons  were  prepared  by 

It  has  now  begun  to  appear  that  the  magnetism  of  arc- melting  weighed  mixtures  of  the  two  metals 

other  Ni-based  alloy  systems,  not  only  of  Ni-Cu,  (both  specified  as  99.  999%  pure)  in  a helium  atmo- 

can  be  characterized  near  cCPlt  in  terms  of  polariza-  sphere  and  were  cold  worked  by  hammering.  The 

tion  clouds  (often  referred  to  as  "magnetic  clus-  samples  machined  from  each  button  were  cylinders 

ters”).  Among  the  evidence  accumulating  to  this  (2-5  mm  diam,  5 mm  long)  whose  edges  were  bev- 

effect,  probably  the  most  convincing  have  again  e(ec*  down  in  order  to  approximate  ellipsoids  (for 

been  the  results  of  elastic  neutron  scattering  exper-  reasonably  uniform  demagnetization).  After  their 
iments  on  weakly  ferromagnetic  alloys.  These  re-  surfaces  were  acid  etched,  the  samples  were  sealed 
suits  have  been  reported  for  Ni-Pd, 24  Ni-Cr, 25  and  in  helium-filled  quartz  tubes,  annealed  for  3 days 

Ni-Rh,2*  and  in  each  case  a marked  upswing  of  the  at  ^00  C,  then  quenched  into  water.  The  nominal 

diffuse  magnetic  cross  section  at  very  small  scat-  compositions  of  the  alloy  samples  are  indicated  in 
tering  vectors  has  signified  the  existence  of  mo-  Table  I.  As  later  inferred  from  our  susceptibility 

ment-density  distributions  extending  over  many  data,  the  ferromagnetic  (presumably  Fe)  impurity 

atomic  volumes,  each  distribution  constituting  a level  in  the  samples  was  probably  well  under  50 

magnetic  cluster.  The  case  of  Ni-Rh  is  particular-  ppm. 

ly  interesting  since  Bucher  et  al . 27  have  claimed,  Magnetization  measurements  on  these  samples 

on  the  basis  of  a specific-heat  and  susceptibility  were  carried  out  with  a temperature-controlled 

study  of  this  system,  that  the  anomalies  observed  vibrating- sample  magnetometer  system  (Princeton 

near  (~63-at.%  Ni)  arise  from  a critical  ex-  Applied  Research  Corp. , Model  No.  150-A).  The 

change  enhancement  accompanied  by  large  spin  Nb-Ti  superconducting  solenoid  (Westinghouse 

fluctuations.  Their  low-temperature  specific-heat  Corp.,  Model  No.  H-5794)  in  this  system  was 

results  were  later  reanalyzed  by  Hahn  and  Wohl-  typically  operated  at  currents  up  to  50  A corre- 

farth23  and  shown  to  be  equally  consistent  with  the  sponding  to  a maximum  field  at  the  sample  of  56 

existence  of  superpara  magnetic  clusters,  for  which  kOe.  The  current-field  specifications  of  the  magnet 

there  also  was  some  qualitative  susceptibility  evi-  were  tested  by  means  of  a calibrated  field- sensi- 

dence  in  alloys  near  Cjj,,.2*  More  recently,  Trip-  tive  Bi-film  resistor  (American  Aerospace  Con- 

lett  and  Phillips30  found  that  the  anomalous  low-  trols  Corp.,  Model  No.  MRA-12);  the  field  homo- 
temperature component  of  the  heat  capacity  of  the  geneity  over  the  sample  length  was  found  to  be  with- 

alloy  NiwRhJS  is  largely  suppressed  (probably  shifted  in  1%.  The  sample-zone  temperature  was  mom- 
to  higher  temperatures)  by  the  application  of  a 38-  tored  by  a GaAs-diode  thermometer  (Lake  Shore 

kOe  field.  They  concluded  that  this  observed  effect  Cryotronics,  Inc. ),  which  we  calibrated  in  zero 

is  in  good  quantitative  agreement  with  the  predicted  field  against  a precalibrated  Ge  resistor  (Cryocal, 
behavior  of  superparamagnetic  clusters  and,  Inc.,  Model  No.  CR-1000)  and  a precalibrated  Pt 

furthermore,  that  it  would  have  required  much  resistor  (Minco  Products,  Inc.,  Model  No.  51061-2) 

larger  fields  if  the  anomalous  heat  capacity  contri-  over  the  ranges  (4.2-40)  K and  (40-260)  K,  respec- 

bution  were  associated  with  spin  fluctuations.  lively.  Each  temperature  calibration  point  corre- 

The  existence  of  magnetic  clusters  in  Ni-Rh  sponded  to  a recorded  setting  of  the  temperature- 

clearly  calls  for  more  complete  documentation  in  control  unit  that  governed  the  current  to  an  elec- 

order  to  establish  whether  or  not  it  is  a significant  trie  heat  exchanger  (temperature  monitored  by 

intrinsic  phenomenon.  In  view  of  this  experimental  another  GaAs  diode),  through  which  vapor  from  the 

need  and  the  continued  interest  in  Ni-Rh,  we  have  liquid-He  reservoir  flowed  into  the  sample  chamber 

carried  out  a systematic  magnetization  study  of  this  at  a controlled  rate.  Since  the  GaAs  diodes  were 

alloy  system  near  the  critical  composition  for  not  entirely  field  insensitive,  calibration  correc- 

ferromagnetism  (cer,|).  Our  preliminary  results  tions  were  determined  as  a function  of  field  at 

for  a few  paramagnetic  Ni-Rh  alloys,  recently  re-  various  fixed  temperatures  set  by  means  of  a He- 

ported  in  brief, 11  appear  to  indicate  that  magnetic  gas  thermometer  bulb  inserted  in  the  sample  cham- 

clusters  do  indeed  play  an  intrinsic  role,  qualita-  ber.  The  maximum  uncertainty  in  the  sample  tem- 

tively  akin  to  their  role  in  Ni-Cu.  In  the  present  perature  over  the  ranges  of  operation  [(4. 2-250)° K, 

paper,  our  complete  results  for  a much  broader  0—56  kOe]  was  about  0.  1 °K;  the  uncertainty  was 

range  of  Ni-Rh  compositions,  extending  through  much  less  at  4.2  °K,  when  the  sample  chamber  was 

ceru,  are  described  and  discussed  in  full  detail.  generally  flooded  with  liquid  He.  The  magnetic- 

pertinent  aspects  of  the  magnetic  critical-point  moment  calibration  of  the  magnetometer  output 

behavior  of  the  weakly  ferromagnetic  alloy  voltage  was  tested  at  different  sensitivity  levels 

Ni,5Rh15,  which  we  have  reported  separately, 32  against  the  saturation  moments  of  various  sized 

are  also  included  in  the  present  discussion.  samples  of  pure  N1  and  Fe  at  4.  2 °K. 33  The  mag- 
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TAHLK  I.  Magnet  li  properties  of  Ni-li!i  alloys. 


Alloy 
(at.T  Ni) 

<r,(0> 

(emu / g) 

<r,(0) 

(emu/g) 

Tc 

CM 

n 

no 

X'  W)* 

r * 

'-cw 

iMB) 

F* 

C) 

CT) 

65.0 

7.40 

11.70 

44 

62 

0.77 

5. 15 

19.7 

0.790 

0.370 

64.4 

5.00 

10.20 

27 

43 

0.91 

4.59 

20.  i 

0.076 

0.328 

64.0 

3.95 

9. 10 

19 

34 

1.05 

4.16 

20.4 

0,  505 

0.  302 

63.2 

1.  70 

7.40 

5 

11 

1.43 

3.76 

22.  7 

0.  43S 

0.257 

62.6 

~0 

5.  80 

-0 

4 

2.00 

2.78 

21.4 

0.  365 

0.227 

62.0 

• • • 

4.95 

... 

_ o 

2 22 

2.29 

20.7 

0. 324 

0.200 

60.0 

• • • 

3.40 

... 

-31 

2.00 

1.61 

21.  1 

0.220 

0. 131 

58.0 

... 

1.55 

... 

-43 

1.76 

0.  826 

23.8 

0.  090 

0.064 

56.0 

... 

0.  96 

... 

-47 

1.63 

0.  501 

23.3 

0.  058 

0.  053 

53.0 

... 

0.52 

-43 

1.33 

0.263 

22.6 

0.033 

0.026 

51.0 

... 

0.  33 

... 

-50 

1.04 

0. 160 

21.7 

0.  022 

0.016 

estimated  uncertainties  (t) 

0.1 

2% 

5% 

0.5 

o 

5rr 

nrr 

3 o 

8T 

13T 

... 

*X'(0)  and  Ccw  in  units  of  lO"’  emu/Oeg  and  UT‘  emu  K/Oeg.  respectively. 


netization  values  for  these  samples  checked  within 
lco,  which  represents  the  absolute  error  in  all  our 
magnetization  measurements.  The  relative  error 
in  the  magnetizations  measured  at  different  tem- 
peratures and  fields  for  a given  mounting  of  a sam- 
ple was  about  an  order  of  magnitude  smaller. 

EXPERIMENTAL  RESULTS  AND  DISCUSSION 

The  magnetization  (cr)  of  each  Ni-Rh  alloy  sample 
was  measured  as  a function  of  magnetic  field  (.//, 
as  corrected  for  demagnetization)  at  temperatures 
(T)  from  4.2  ’K  up  to  about  250  °K.  Our  results 
for  a vs  H at  4.  2 °K  are  illustrated  in  Fig.  1.  The 
curves  for  the  alloys  of  63.2-at.^  or  more  Ni  show 
a spontaneous  magnetization  ta0)  clearly  indicative 
of  ferromagnetism.  The  Curie  points  (Tc)  of  these 
alloys  were  determined  from  isothermal  a2  vs  H/o 
plots  of  the  data;  their  values,  together  with  those 
of  cr0  extrapolated  to  0 3K,  are  listed  in  Table  I.  A 
similar  plot  for  the  62.6-at.  % Ni  sample  gives  Tc 
=<0’K,  indicating  that  this  alloy  is  essentially  at  the 
critical  composition  for  ferromagnetism  (c=rlt ). 
Moreover,  the  cr2  vs  H/o  plots  near  Tc,  which  are 
extremely  sensitive  to  any  macroscopic  inhomo- 
geneities, revealed  no  metallurigical  complications 
in  our  alloy  samples. 

For  each  alloy  over  its  paramagnetic  range  of 
temperatures,  the  initial  susceptibility  lx0)  was  de- 
termined from  da/dH  extrapolated  to  H-  0;  its  re- 
ciprocal values  (yj1)  are  shown  plotted  versus  T in 
Fig.  2.  In  their  concave-downward  shape  Iparticu- 
larly  evident  at  low  temperatures),  these  curves 
agree  with  the  earlier  Nl-Rh  results  of  Cottet  et 
nl  ,29  They  differ  qualitatively  from  the  flat  concave- 
upward  curves  normally  predicted  from  simple 
itinerant-electron  band  theory.  * In  this  sense,  the 
Xo'  vs  T curves  for  Ni-Rh  resemble  those  for  para- 
magnetic Ni-Cu  near  where  the  latter  were 


each  found  to  be  analyzable  into  a Curie- Weiss 
component  and  a weak  temperature-independent 
component,  attributed  respectively  to  superpara- 
magnetic  clusters  and  itinerant-electron  band 
polarization.  a 

Following  an  analogous  procedure  for  Ni-Rh  al- 
loys near  ccrlt,  we  assume  that  the  measured  mag- 
netization may  be  expressed  as 

<j(H,T)  = <j'1W,T)  + Hx\T)  , (1) 

where  eta  represents  the  magnetization  arising  from 
the  spin  alignment  of  clusters  and  \ is  a band- 
polarization  susceptibility,  which  is  expected  to  be 
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FIG.  1.  Magnetization  vs  field  -at  4.2*K  for  Ni-Rh 
alloys  (compositions  in  at.  % Ni). 
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FIG.  2.  Inverse  initial  susceptibility  vs  temperature 
for  Ni-ilh  alloys  (compositions  in  at.  % Ni). 


only  mildly  temperature  dependent  and  negligibly 
dependent  on  laboratory  fields.  Consequently,  the 
spontaneous  magnetization  of  a ferromagnetic  al- 
loy below  Tc  will  derive  entirely  from  the  zero- 
field  alignment  of  cluster  moments,  i.  e. , o0(T) 

= acll0,  T),  whereas  the  initial  susceptibility  above 
Tc  land  at  all  temperatures  for  a paramagnetic  al- 
loy) will  include  a band  polarization  term  as  well 
as  a superparamagnetic  cluster  term,  the  latter 
of  which  is  taken  to  be  Curie- Weiss-like,  i.e., 

xoin-x^n+x'w) 

= cc„Ar-e)+x\r)  . (2) 

In  the  high-field  limit,  achieved  by  extrapolation 
from  laboratory  fields,  it  follows  that 

T)  = lim[cr(tf , T)  - 13) 

#-» 

and 

X.\T)  = lim(d<j/dH)T  , 14) 

«»» 

where  a,  represents  the  saturation  limit  of  the  clus- 
ter spin  alignment. 

In  order  to  separate  our  y0(r)  data  for  Ni-Rh  into 
the  components  indicated  in  Eq.  (2),  we  took  as  an 
approximate  low-temperature  value  for  the  dif- 
ferential susceptibility  at  4.2  °K  extrapolated  to  in- 
finite field,  as  suggested  by  Eq.  (4).  The  value  thus 
estimated  was  ~2x  10*5  emu/Oeg  for  all  the  para- 
magnetic alloys.  This  starting  value  for  x’  was  then 
adjusted  for  each  of  these  alloys  until  Ixo-x’)*' 


vs  T became  linear  over  a range  of  low  tempera- 
tures (i.e.,  up  lo~20t,K).  In  order  that  this  lin- 
earity continue  up  to  moderate  temperatures  (~  100 
°K),  it  was  found  necessary  in  each  case  to  let  x* 
have  a weak  temperature  dependence  of  the  form 

x'(^)  = [x’*‘(0)  + arz]-1  (5) 

with  a positive.  At  still  higher  temperatures, 
where  the  latter  procedure  began  to  overestimate 
the  variation  in  x . we  simply  extended  the  linear 
dependence  of  (Xo-x')"1  on  T and  computed  back- 
wards to  obtain  the  \ ' values  corresponding  to  the 
measured  values  of  x0-  Our  results  for  the  tem- 
perature dependence  of  x*i'  = (Xo  ~ xV'  and  of  x'1 
are  plotted  in  Figs.  3(a)  and  (b),  respectively;  the 
values  of  x'(0)  are  listed  in  Table  I.  For  the  fer- 
romagnetic alloys,  since  the  above  analysis  can- 
not be  used  below  Tc,  we  adjusted  x (0)  to  the  val- 
ues listed  in  Table  I and  assumed  that  x *l(T)  is  dis- 
placed upward  uniformly  from  the  curve  determined 
for  the  critical-composition  (62.  6-at.  % Ni)  alloy, 
as  illustrated  in  Fig.  3(b).  The  linear  plots  thus 
obtained  for  xIi‘(T)  above  Tc,  which  are  presented 
in  Fig.  3(a),  are  rather  insensitive  to  the  exact 
shape  assumed  for  owing  to  the  fact  that 

is  only  a small  part  of  xo  over  this  temperature 
range. 

For  each  Ni-Rh  alloy,  the  values  of  Ccw  and  9 
obtained  from  the  Curie- Weiss  behavior  of  x^lT), 
as  expressed  in  Eq.  (2),  are  listed  in  Table  I.  Al- 
lowing that  Xd  is  associated  with  superparamag- 
netic  clusters,  we  can  write  for  the  Curie-Weiss 
constant 

CCw  = .Yc*(i i*J>/3fe  , (6) 

where  S is  the  number  of  atoms  per  gram,  c*  the 
cluster  concentration  (number  per  atom),  (pmi)  the 
mean  squared  value  of  the  magnetic  moment  per 
cluster,  and  k the  Boltzmann  constant.  We  further 
allow  that  the  saturation  magnetization  o,(0)  ob- 
tained by  high- field  extrapolation  of  o(/f,  0)  - 
as  indicated  in  Eq.  (3),  represents  tha  parallel  spin 
alignment  of  all  the  magnetic  clusters  in  the  alloy 
and  can  therefore  be  expressed  as 

a,(0  )»Nc*<P*>.  (7) 

Typical  high-field  extrapolations  for  the  Ni-Rh  al- 
loys, involving  our  low- temperature  magnetiza- 
• tion  data  (extended  trivially  to  0 °K)  and  the  values 
deduced  for  x’lO),  are  illustrated  in  Fig.  4.  The 
values  thus  determined  for  o,(0)  are  listed  in  Table 
I.  From  a combination  of  Eqs.  (6)  and  (7),  we  ob- 
tain the  “average”  quantities  defined  as 

P*  a(p*l)/(g*>  * 3/.-Ccw/o,(0)  , (8a) 

<T*  *c*(p*)/7l*  = oJ(0)/3.VfeCcw  . (8b) 

(Note  that  a distribution  in  cluster  moment  size  will 
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FIG.  3.  (al  Inverse  ol  Xei5Xo~X'  and  (b>  inverse  of 
X',  as  functions  of  temperature,  for  Ni-Rh  alloys  (com- 
positions in  at. % Ni). 


affect  the  meaning  of  p*  and  cm,  rendering  them 
respectively  higher  and  lower  than  the  true 
average  values,  (p*>  and  c“.)  Substituting  our  ex- 
perimentally derived  values  for  Ccw  and  0,(0)  into 
these  expressions,  we  calculated  the  m*  and  c* 
values  listed  in  Table  I.  From  these  results  we 
see  that  the  cluster  moment  has  enormous  val- 
ues [-  (20-24)pa]  that  are  remarkably  constant  over 
the  composition  range  of  study.  Contrastingly,  the 
cluster  concentrations  c*  are  very  small  and  de- 


n 

crease  rapidly  with  decreasing  at.'l  Ni,  passing 
smoothly  through  the  critical  composition.  Qualita- 
tively, these  results  for  the  average  moment  and 
concentration  of  magnetic  clusters  in  Ni-Rh  are 
very  similar  to  'hose  previously  obtained  for12 
Ni-Cu  and34  Ni-V  by  the  same  method  of  data 
analysis,  which  yielded  giant  g*  values  of  - I0p.a 
and-40jia,  respectively,  near  ecrlt. 

Evidence  for  giant  magnetic  clusters  in  the  weak- 
ly ferromagnetic  Ni65Rh35  alloy  was  also  extracted 
from  detailed  ot H,  T)  data  in  the  vicinity  of  Tc,  as 
we  have  already  reported  separately.32  From  these 
data,  the  spontaneous  magnetization  just  below  Tc 
and  the  initial  susceptibility  just  above  Tc  were  de- 
termined and  their  temperature  dependences  de- 
scribed respectively  as  o0t,T)=  iual  1 - T/Tcl‘  and 
\0(T)  = (wi0A0)l  1 - T/Tc i~r,  where  |3  and  > are  crit- 
ical exponents,  and  w0  and  inQ/fi0  may  be  regarded 
as  critical  coefficients.  Our  results  for  the  latter 
were  compared  with  the  theoretical  predictions 
that  />t0/a0(O)  = it  | and  nji0/kTc  = n2,  in  which  the  ele- 
mentary moment  is  normally  taken  to  be  the 
average  moment  per  atom  and  the  numbers  it,  and 
it2  are  both  slightly  larger  than  unity  for  various 
classic  models  (mean  field,  Ising).  While  our  ex- 
perimental value  for  it,  obeyed  these  predictions 
quite  well,  the  value  for  »2  was  only  -0.00"  when 
Mq  was  equated  to  0oiO)/.V=O.  lga.  Indeed,  agree- 
ment with  the  near-unity  values  predicted  for  n2  re- 
quired that-we  set  g0  = 22 (ia,  which  is  essentially 
the  p*  value  for  this  alloy  given  in  Table  I.  Thus, 
it  was  concluded  that  the  elementary  dynamical  mo- 
ments involved  in  the  Curie-point  transition  are 
those  of  giant  magnetic  clusters  rather  than  of  in- 


FIG.  4.  Variation  of  <i-\'H  <for  T— 0*K)  with  W'  for 
Ni-Rh  alloys  (compositions  in  at.  T N ). 
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dividual  atoms.  We  further  showed  that  the  sami 
conclusion  can  be  reached  analogously  from  the 
C iric-point  properties  of  dilute  Pd Fe,  the  classic 
g at- moment  alloys— whose  critical  exponent  val- 
ues, incidentally,  also  are  anomalous  and  similar 
to  those  of  the  Ni-Rh  alloy. 

As  was  done  in  the  case  of  Ni-Cu, 12  let  us  exa- 
mine the  possibility  that  the  magnetic  cluster  con- 
centrations in  Ni-Rh  correlate  with  the  statistical 
occurrences  of  a particular  local  chemical  compo- 
sition. For  simplicity,  consider  the  probability  P„ 
that  any  lattice  site  of  a fee  NicRh,^  alloy  (assumed 
to  be  a random  solid  solution)  is  occupied  by  a Ni 
atom  whose  12  nearest  neighbors  include  n or  more 
Ni  atoms.  In  the  extreme  case  of  n = 12,  we  have 
P^-c11,  whose  values  for  the  alloys  studied  are 
listed  in  Table  I.  The  variation  of  Pl2  with  alloy 
composition  is  clearly  quite  similar  to  that  of?*, 
although  at  the  Ni-rich  end  of  this  range  the  values 
of?*  appear  to  be  rising  towards  Pn  = c15 
♦ I2e12(l  -c),  whose  values  exceed  those  of  Pl2  by 
nearly  an  order  of  magnitude.  Despite  the  over- 
simplicity of  these  considerations,  they  are  prob- 
ably correct  in  indicating  that  the  giant  magnetic 
clusters  in  Ni-Rh  are  nucleated  in  local  regions 
that  are  extremely  Ni  rich.  In  the  case  of  Ni-Cu, 

?*  was  originally  found12  to  lie  between  Pl0  and  Pn, 
but  it  was  recently  shown20  in  a more  rigorous  sta- 
tistical treatment  of  the  atomic  clustering  that?* 
follows  closely  the  values  of  Pa.  Thus,  the  local 
conditions  for  magnetic  cluster  nucleation  in  Ni-Rh 
and  Ni-Cu  appear  to  be  fairly  alike.  In  either  case, 
according  to  the  neutron  diffraction  work  on  Ni- 
Rh,  28  and  Ni-Cu.  11  the  average  moment  density  in- 
side a magnetic  cluster  decreases  gradually  over 
a radial  distance  of  many  lattice  spacings  and  thus 
extends  considerably  outside  the  central  core  of 
nearest-neighbor  Ni  atoms. 

Whatever  may  be  the  underlying  mechanisms  for 
magnetic  cluster  formation  in  Ni-Rh  (which  will  be 
commented  upon  later),  it  is  clear  from  our  results 
that  giant  magnetic  clusters  persist,  as  superpara- 
magnetic  entities  in  extremely  dilute  concentra- 
tions, well  into  the  paramagnetic  composition 
range.  The  negative  values  of  the  Curie- Weiss 
temperature  (0)  for  the  weakly  paramagnetic  Ni- 
Rh  alloys,  similar  to  (though  larger  in  magnitude 
than)  the  negative  0 values  for  paramagnetic  Ni- 
Cu,  12  and  Ni-V,34  may  reflect  local  anisotropy  ef- 
fects38 as  much  as  a predominance  of  antiferromag- 
netic over  ferromagnetic  exchange  interactions, 
presumably  of  the  indirect  Ruderman-Kittel-Kas- 
uya-Yosidatype.19  As  ccrlt  is  approached  from  the 
paramagnetic  side  and  the  cluster  concentration  in 
Ni-Rh  increases,  0 rises  sharply  to  large  positive 
values,  probably  as  the  outcome  of  increased  fer- 
romagnetic (direct  overlap)  interactions  between  ad- 
jacent clusters.  Beyond  cer„,  in  the  weakly  ferro- 
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magnetic  Ni-Rh  alloys,  the  spontaneous  magnetiza- 
tion o0(0)  is  much  smaller  than  the  corresponding 
saturation  value  0,(0),  although  their  ratio  does  ap- 
proach unity  with  increasing  at. 'o  Ni,  as  seen  in 
Table  I.  This  behavior  shows  that  in  zero  external 
field  many  of  the  magnetic  clusters  in  these  alloys 
are  not  ferromagnetically  aligned  even  at  very  low 
temperatures— which  is  consistent  with  the  fact  that 
specific-heat  anomalies  of  the  type  ascribable  to 
superparamagnetic  clusters29  are  observed  in  the 
weakly  ferromagnetic  as  well  as  the  strongly  para- 
magnetic Ni-Rh  alloys.2’ 

According  to  our  interpretation,  the  onset  of  fer- 
romagnetism in  Ni-Rh  arises  from  the  exchange 
coupling  between  magnetic  clusters  and  is  mani- 
fested at  Tc  in  the  divergence  of  \ci>  the  Curie- 
Weiss-like  cluster  component  of  the  initial  sus- 
ceptibility, as  indicated  in  Fig.  3(a).  The  other 
component  of  the  initial  susceptibility  (i.e. , \ ) 
exhibits  only  a mild  T2  deviation  from  its  zero- 
temperature  value  and,  as  further  seen  from  Fig. 
3(b),  the  variation  of  x"1  at  higher  temperatures 
seems  to  approach  a weak  linear  dependence.  A 
very  similar  behavior  was  recently  deduced  for  \ 
vs  T in  Ni-V  alloys  close  to  the  critical  composi- 
tion,34 for  which  x\o)  is  comparable  in  size  to  our 
Ni-Rh  results.  Values  obtained  for  \’  in  Ni-Cu  al- 
loys near  ccrlt  are  about  an  order  of  magnitude 
smaller  and  were  assumed  to  have  negligible  tem- 
perature dependence12;  they  were  attributed,  as 
was  mentioned  earlier,  to  weak  band  polarization 
effects.  In  the  case  of  Ni-Rh,  the  x *1  vs  T curves 
in  Fig.  3(b)  have  exactly  the  shape  predicted  for 
such  effects  when  the  exchange  enhancement  is 
large  and  accompanied  by  spin  fluctuations. 36  An 
analogous  explanation  was  given  for  the  similar 
x'*‘  vs  T curves  for  Ni-V. 37 

Our  x^O)  results  for  Ni-Rh  are  plotted  versus 
alloy  composition  in  Fig.  5,  together  with  dadH 
values  taken  at  various  fields  from  the  a vs  H 
curves  for  4.  2 °K  in  Fig.  1.  The  peaks  in  the  da/ 
dH  curves  at  the  critical  composition  (~63-at.‘o 
Ni),  which  are  most  prominent  at  the  lowest  fields, 
are  clearly  the  vestiges  of  the  divergence  of  da/dH 
at  zero  field.  Hence,  they  relate  more  closely  to 
the  magnetic  cluster  component  of  the  initial  sus- 
ceptibility (i.e.,  Xei)  than  to  ttie  bant*  polarization 
component  (x*),  and  this  appears  to  be  ti;ue  even  at 
the  highest  fields  where,  according  to  Eq.  (4),  da! 
dH  should  be  extrapolating  to  x*.  It  is  evident  from 
Fig.  5 that  the  high-field  approach  of  da/dH  to  x* 
is  extremely  slow  in  all  these  alloys,  not  only  in 
those  very  close  to  This  behavior  quite  pos- 

sibly derives  from  a broad  distribution  in  the  mag- 
netic moments  of  the  individual  clusters,  in  which 
case  the  clusters  with  relatively  small  moments 
continue  to  contribute  to  the  low-temperature  dif- 
ferential susceptibility  up  to  very  high  fields. 
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FIG.  5.  Open  circles:  do/ilH  at  various  fields  (in  kOe) 
at  4.2'K;  closed  circles:  x ' <0) , as  functions  of  Ni-Rh 
alloy  composition. 


Cluster  moment  distributions  in38  Ni-Cu  and34  Ni-V 
have  been  proposed  on  the  basis  of  similar  experi- 
mental evidence. 

The  variation  of  our  results  for  x'lO)  is  shown  in 
Fig.  6 with  respect  to  the  entire  Ni-Rh  composition 
range.  On  this  scale,  the  values  for  x^O)  define 
a sharp  but  finite  peak  centered  near  the  critical 
composition.  Their  rapid  decrease  on  the  Rh-rich 
side  of  the  peak  extends  smoothly  through  the  points 
representing  some  previous  high-field  susceptibility 
data, 37  which  at  these  compositions  can  be  safely 
assumed  to  contain  virtually  no  magnetic  cluster 
contribution.  The  descent  of  x,(0)  on  the  Ni-rich 
side  of  the  peak  appears  to  be  even  faster,  probably 
reaching  very  low  values  for  most  of  the  ferromag- 
netic range,  as  indicated  in  the  figure.  This  be- 
havior suggests  that  the  band  polarization  effects 
that  produce  x',  though  not  responsible  for  the  on- 
set of  ferromagnetism,  become  highly  exchange  en- 
hanced in  the  vicinity  of  c^,,.  It  seems  reasonable 
to  associate  this  strong  exchange  enhancement  with 
local  (presumably  fairly  Ni  rich)  regions  that  are 
nearly  unstable  towards  the  formation  of  local 
(cluster)  moments.  As  implied  by  our  x'  results 
in  Fig.  6,  these  local  regions  grow  and  approach 
critical  enhancement  as  the  alloy  composition  in- 
creases in  at.%  Ni,  but  at  caix  they  suddenly  be- 
come exchange-polarized  by  the  spontaneously 
aligned  moments  of  the  stable  magnetic  clusters 


and,  as  a result,  their  contribution  to  the  sus- 
ceptibility (x')  drops  precipitously.  Spin  fluctua- 
tions can  be  expected  to  occur  in  these  regions  of 
strong  local  enhancement  and  in  fact,  according 
to  a recent  resistivity  study, 39  they  appear  to  pro- 
vide the  dominant  scattering  mechanism  in  Ni-Rh 
alloys  near  ccrll.  In  the  case  of  Ni-V,  where  x’ 
attains  values  comparable  to  those  in  Ni-Rh,  the 
observed  temperature  dependence  of  the  resistivity 
was  also  found  to  be  characteristic  of  spin-fluctua- 
tion scattering34;  as  in  Ni-Rh,  any  scattering  from 
the  superparamagnetic  clusters  seems  to  be  com- 
pletely obscured.  However,  in  the  case  of  Ni-Cu, 
where  x'  is  very  much  weaker,  the  measured  re- 
sistivity versus  temperature  curves  indicate  that 
the  predominant  scattering  is  from  superparamag- 
netic clusters,40  as  was  recently  elucidated  theo- 
retically.21 

Having  attributed  the  x(0)  peak  in  Fig.  6 to  re- 
gions of  near-critical  exchange  enhancement,  we 
would  be  consistent  in  regarding  each  statistical 
Ni-rich  local  region  that  constitutes  the  core  of  a 
magnetic  cluster  in  Ni-Rh  as  having  reached  and 
exceeded  critical  enhancement.  There  would  con- 
sequently exist  stable  local  moments  on  all  the  Ni 
atoms  (and  perhaps  the  Rh  atoms  as  well)41  within 
each  cluster  core.  From  this  point  of  departure  in 
the  case  of  Ni-Cu,  Roth  has  proceeded  to  show  ana- 
lytically that  the  Ni  atoms  outside  the  core  of  a 
magnetic  cluster  will  become  magnetically  polar- 
ized but  to  a decreasing  extent  with  increasing  dis- 
tance from  the  core.  17  This  picture  of  a giant  mag- 
netic cluster  or  polarization  cloud  in  Ni-Cu  agrees 
with  the  spin  density  distribution  deduced  from 
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FIG.  6.  Closed  circles:  x ' (0)  vs  Ni-Rh  alloy  compo- 
sition; open  circles:  high-field  susceptibility  data  from 
Ref.  27. 
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neutron  diffraction  data  and  has  subsequently  been 
described  by  Garland  and  Gonis  in  a sclf-consisti  lit 
theoretical  treatment  of  the  entire  magnetic  cluster 
(core  and  all).20  The  cooperative  magnetic  pro- 
cesses invoked  in  the  latter  work  are  presumably 
also  responsible  for  the  giant  magnetic  clusters  in 
Ni-Rh,  but  this  remains  to  be  investigated. 
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Giant  Moments  in  Paramagnetic  Pd-Ni  Alloys* 

D.  SAIN  and  J.  S.  KOUVEL,  Univ.  of  Illinois,  Chicago  — 
For  Pd-Ni  alloys  approaching  the  critical  composition  for 
ferromagnetism  (^2.5  at./!  Ni)  the  initial  susceptibility 
as  a detailed  function  of  temperature  T can  be  resolved 
into  a Curie-Weiss  component,  Xcv  = Ccv/(T-0),  and  a com- 
ponent x'  vhose  weak  dependence  on  T is  characteristic  of 
unstable  (fluctuating)  moments.  Up  to  'VL.8  a.t.%  Hi, 
where  x'(0°K)  reaches  a maximum,  CQ  increases  slowly  and 
is  consistent  with  giant  moments  associated  with  the  iron 
impurity  atoms  (a<80  ppm)  in  the  alloys.  Above  1.8  at.)! 
Ni,  both  Ccw  and  the  saturation  magnetisation  at  b.2°K 
rise  very  rapidly  and  are  attributable  to  an  increasing 
concentration  of  giant-moment  clusters  centered  at  small 
statistical  groups  of  neighboring  Ni  atoms.  The  values 
of  0 are  extremely  small  (<0.1°K)  up  to  2.1  at./!  Ni  and 
become  increasingly  positive  as  the  critical  composition 
is  approached  more  closely. 
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Detailed  magnetization-field-temperature  measurements  on  C.'-rich  Fe-Cr  alloys  near  their  ferromagnetic  Curie  points 
reveal  anomalous  values  for  the  critical  exponents.  From  the  values  for  the  critical  coefficients,  it  is  deduced  that  the 
alloys  are  giant-moment  systems,  similar  to  Pd-Fe  and  Ni-Rh 


I.  Introduction 

The  bulk  magnetic  properties  [1.2]  of  Fe-Cr 
alloys  containing  greater  than  30at.%Fe  indi- 
cate that  the  alloys  are  normal  ferromagnets, 
and  extrapolation  of  the  data  to  higher 
chromium  concentrations  suggests  that  long- 
range  ferromagnetic  order  should  disappear  in 
the  region  15-20  at.%  Fe.  Neutron-diffraction 
experiments  [3]  imply  that  in  this  critical 
concentration  region  the  magnetic  moments  will 
be  confined  primarily  to  the  iron  atoms,  with  a 
magnitude  that  will  vary  profoundly  as  a func- 
tion of  the  local  atomic  environment.  By 
analogy  with  the  Ni-Cu  system  [4],  we  would 
thus  anticipate  the  presence  of  giant  moments  in 
the  critical  concentration  range  centered  about 
iron-rich  clusters. 

Earlier  studies  [5-8]  of  the  magnetic  proper- 
ties of  Fe-Cr  alloys  near  the  critical  concen- 
tration do  indicate  that  the  distribution  of  mag- 
netization is  inhomogeneous.  However,  it  is  not 
clear  whether  the  various  results  are  self  con- 
sistent and  whether  the  alloys  can  be  described 
as  ferromagnetic  or  mictomagnetic.  The 
presence  of  short-range  chemical  inhomogenei- 
ties in  the  Fe-Cr  system  [3]  complicates  this 
problem,  inasmuch  as  no  two  nominally-iden- 
tical  alloys  prepared  on  separate  occasions  will 
ever  have  precisely  the  same  atomic  distribution 
and,  therefore,  the  same  magnetic  properties. 

In  the  present  investigation,  we  have  deter- 
mined the  detailed  magnetic  properties  of  Fe-Cr 
alloys  containing  30.  25,  and  20at.fTFe  in  the 

• Work  supported  by  the  US  Energy  Research  and  De- 
velopment Administration 
♦Supported  in  part  by  NSF  and  ONR 


vicinity  of  the  Curie  temperature  Tc.  The  pri- 
mary aim,  following  the  study  of  Kouvel  and 
Comly  [9],  has  been  to  determine  the  magnetic 
entity  taking  part  in  the  ordering  process  and 
relate  this  to  the  anticipated  ‘giant  moment’. 

2.  Experimental  procedure  and  results 

The  alloys  were  prepared  by  arc  melting  the 
requisite  amounts  of  99.99  wt.%  pure  iron  and 
chronium  in  a helium-argon  atmosphere.  After 
suitable  heat  treatment  [2],  samples  were  ground 
in  the  form  of  prolate  ellipsoids  (— 5mmx 
— 2 mm)  and  reannealed  at  900°C  for  one  day. 
Chemical  analysis,  and  X-ray  lattice  parameter 
and  density  determinations  all  indicated  that, 
within  error,  the  actual  and  nominal  compositions 
were  the  same. 

A series  of  magnetic  isc.herms  [specific 
magnetization  <r(emu/g)  as  a function  of  inter- 
nal field  H(Oel]  was  determined  at  l K.  intervals 
by  the  Faraday  method  in  an  apparatus  already 
described  [10],  As  a function  of  both  increasing 
and  decreasing  applied  field  in  the  range  0.25- 
I4k0e,  the  data  yielded  a set  of  smooth  curves 
when  plotted  as  cr:  versus  Hlcr. 

Extrapolation  of  these  curves  to  the  axes 
gave  the  inverse  initial  susceptibility  Yo~'(T  > Tc) 
and  the  square  of  the  spontaneous  magnetiza- 
tion cr's(T  < Tc).  Analysis  of  the  temperature 
dependence  of  these  quantities  determined  Tc 
and  the  critical  exponents  y and  /?  defined  by: 

«i=m  o‘Va  and  = (ho/m„)e\  (1.2) 

where  g = |l  - 77 Ts |,  and  m0J  and  h0lm„  are 
critical  coefficients.  The  critical  isotherm 
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Table  I 

Magnetic  critical-point  parameters  for  the  alloys  listed 


Fe  «,Cr.« 

Fe,.Cr„ 

Fe„Cr«, 

Nu.Rh,. 

Pd-.Fe,. 

Ni 

T.(K) 

2501 1 )• 

154(1) 

64(2) 

44  2 

50.2 

627  4 

tr„  (emu/g) 

57.6 

42 

18 

7 40 

6 85 

58  57 

0 548 

0 40 

0.17 

0.098 

0 130 

0.616 

0 

0.444 1 ) 

047(3) 

0.52(1) 

0.476 

0 464 

0 378 

y 

1.82(1) 

2.03(4) 

2.12(1) 

1.50 

1.40 

1.34 

8 

5.15(3) 

5.32(3) 

5.11(3) 

4.15 

4.06 

4.58 

5.18 

5.30 

5 10 

4.15 

4.02 

4.54 

m„(e  mu/g) 

63.6 

47.2 

26.8 

8.66 

8.25 

83.3 

6„(kOe) 

509 

236 

28 

48.3 

132 

15720 

mj 

1 10 

1.12 

1.3 

1.17 

1 20 

1.42 

liohj  k Tt 

0.075 

0044 

0 0041 

0.0072 

0.023 

1.04 

12 

16 

55 

22 

9 

- 

c*  (%) 

4.6 

2.6 

0.31 

0.45 

1.4 

~ 

'Numbers  in  parentheses  represent  estimated  uncertainties  in  least  significant  figure 


(T  = TJ  was  then  plotted  as  log  < t versus  log  H 
to  give  the  exponent  6(rr  * H1*). 

The  values  obtained  for  the  various 
parameters  are  listed  in  table  1.  together  with 
comparative  results  for  Ni  [II],  a Pd^Fe,., 
alloy  [9|,  and  a Ni^,Rh„  alloy  [12].  The  unusual 
exponent  values  found  in  this  study,  particularly 
for  y.  may  be  related,  by  analogy  with  previous 
results  |9.  12],  to  an  inhomogeneous  distribution 
of  magnetization.  In  all  cases,  the  exponents 
conform  to  the  scaling  relation  5=1+  yip,  as 


I 


Fig  I Plot  of  some  normalized  isotherms  (see  lest)  near  T, 
for  an  Fe„Cr».  sample. 


shown  by  the  value  of  5 calculated  in  this 
manner  from  the  y and  (3  values  given  in  table  I. 
The  present  data  also  obey  the  homogeneous 
equation  of  state,  as  shown  by  fig.  1 where  some 
of  the  results  for  the  Fe.,Cr-<  sample  are  plotted 
[9.12]  as  versus  hlm[  = (Hlcr)fy]. 

3.  Discussion 

The  intercepts  of  the  curves  in  fig.  I with  the 
axes  determine  the  coefficients  rn;,  and  hjm„  in 
eqs.  (I)  and  (2),  and  these  may  be  combined  to 
yield  m0  and  h„  which  are  listed  in  table  I.  Also 
listed  are  the  normalized  quantities  mjcr„  and 
p„hJkTc,  where  <rn  is  the  spontaneous  mag- 
netization at  absolute  zero  and  /!„  is  the  cor- 
responding mean  moment  per  atom.  The  uncer- 
tainty in  the  <r„  values  for  the  20  and  25at.'T-Cr 
samples  arises  from  the  pronounced  field 
dependence  of  the  magnetization  even  at  low 
temperatures  which  made  a reliable  ex- 
trapolation to  H = 0 difficult.  Although  the 
values  of  are  close  to  previous  results  and 

theoretical  estimates  [9.  12],  the  values  of 
p.„hJkTc  are  more  than  an  order  of  magnitude 
less  than  theoretical  values,  as  in  the  case  of 
Pdg,ftFe,4  [9]  and  Ni„,Rh„  [12].  Following  pre- 
vious work  (9.  12],  we  define  /x*  as  the  effective 
elementary  moment  taking  part  in  the  ordering 
process,  and  equate  n*hJkT<  with  the  theoreti- 
cal estimate  (—  1.6  [12])  to  obtain  values  of  m* 
the  giant  moment  which  are  given  in  table  1. 
Furthermore,  if  the  concentration  of  these  en- 
tities is  c*.  then  ii„=  n*c*  and  values  of  c* 
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obtained  in  this  way  arc  also  listed  in  table  1. 
Itecause  of  uncertainties  in  the  extrapolation 
procedures  and  the  theoretical  estimate  of 
n*hJkT c.  the  values  of  n*  and  c*  should  be 
considered  approximate.  However,  a general 
trend  appears  in  that,  as  -the  critical  concen- 
tration is  approached  from  the  ferromagnetic 
side,  the  giant-moment  clusters  become  larger  in 
magnitude  and  smaller  in  concentration. 

From  these  results,  we  conclude  that  the 
Curie-point  transition  in  weakly  ferromagnetic 
Fe-Cr  alloys  involves  the  magnetic  order- 
disorder  of  a dilute  concentration  of  giant- 
moment  clusters.  Moreover,  the  critical  ex- 
ponents that  characterize  this  transition  have 
anomalous  values  similar  to  those  of  other 
giant-moment  alloy  systems. 
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The  isothermal  magnetization  processes  in  the  alloy  Cu81Mni,  are  found  to 
undergo  a qualitative  change,  involving  the  sudden  appearance  of  magnetic 
hysteresis,  at  the  temperature  (~  85°K.)  at  which  the  initial  susceptibility 
has  a sharp  maximum. 


THE  COMPLEX  magnetic  properties  of  Cu-Mn  and 
certain  other  transition-group  alloy  systems  have  been 
attributed  to  a correspondingly  complex  magnetic 
state  known  variously  as  a mictomagnetic  or  spin- 
glass  state.1  Whereas  the  mictomagnetic  picture  for 
these  atomically  disordered  alloys  stresses  an  inhomo- 
geneous coexistence  of  ferromagnetic  and  antiferro- 
magnetic local  spin  alignments,  the  spin-glass  picture 
emphasizes  a nearly  random  distribution  of  atomic 
spin  orientations.  At  low  temperatures,  both  pictures 
allow  for  strong  local  anisotropy  forces  which  “freeze 
in”  a static  spin  configuration  with  no  long-range  mag- 
netic order. 

The  important  role  in  Cu-Mn  of  some  type  of 
local  magnetic  anisotropy  was  first  seen  in  the  un- 
usual displaced  hysteresis  loops  which  result  from 
cooling  the  alloys  in  a magnetic  field.2  Furthermore, 
when  an  alloy  in  this  metastable  Held-cooled  state  is 
warmed  up,  there  is  an  enormous  increase  in  mag- 
netic hysteresis  losses  as  the  measured  loops  gradually 
lose  their  asymmetry.2  In  apparent  contrast  with  this 
gradual  process,  the  initial  flow-field)  susceptibility  of 
Cu-Mn  was  recently  found  to  exhibit  a cusp-like 
maximum.3  For  a given  alloy  composition,  this  sharp 
phenomenon  occurs  at  essentially  the  temperature  at 
which  the  effects  produced  by  field-cooling  disap- 
pear.2,3  Although  the  sharp  susceptibility  peak  may 

* This  work  was  supported  in  part  by  the  National 
Science  Foundation  and  the  Office  of  Naval  Re- 
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suggest  a Neel-point  transition  to  an  antiferromagnet- 
ically  ordered  state,  the  neutron  diffraction  evidence 
for  any  type  of  long-range  magnetic  order  in  Cu-Mn 
is  quite  negative.4  The  current  interpretation  of  this 
sharp  behavior  is  that  it  marks  the  rapid  freezing  in  of 
a complex  glass-like  spin  configuration  of  the  type 
mentioned  above.  However,  this  interpretation  is  not 
definitive,  and  many  interrelated  properties  need  to 
be  explored  before  the  magnetic  states  of  Cu-Mn  can 
be  properly  understood. 

In  this  paper,  we  report  on  a magnetization 
study  of  a concentrated  Cu-Mn  alloy.  The  isothermal 
magnetization  processes  are  found  to  undergo  a qual- 
itative change,  involving  the  sudden  appearance  of 
magnetic  hysteresis,  at  the  temperature  at  which  the 
initial  susceptibility  has  a sharp  maximum.  Thus,  the 
sharp  susceptibility  peak  is  seen  as  one  of  several 
manifestations  of  an  unusual  magnetic  transition. 

For  this  study,  a piece  was  taken  from  an  alloy 
ingot  (prepared  by  induction-melting  under  argon)  of 
nominal  atomic  composition  Cu^Mn^,  and  was  cold- 
worked  by  hammering.  A cylindrical  sample  was  cut 
from  this  ingot  piece,  annealed  for  three  days  at 
800°C  and  water-quenched,  and  tlien  aged  by  slow- 
cooling  from  350°C  to  room  temperature  at  a rate  of 
lOdeg/day.  The  magnetization  measurements  were 
made  with  a vibrating-sample  magnetometer  at  tem- 
peratures down  to  4.2°K  and  in  fields  up  to  56  kOe. 
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FlG  . I . Magnetization  vs  field  curves  for  CusiMn,g 
cooled  to  various  temperatures  (in  °K)  in  zero  field; 
each  curve  was  measured  isothermallv  with  increas- 
ing field.  Insert:  initial  susceptibility  vs  temperature 
as  determined  from  the  initial  slopes  of  the  various 
isotherms. 


FlG  2.  Magnetization  vs  field  for  Cu9|Mni9  at  4.2°K. 
Solid  curves  measured  after  cooling  in  zero  field,  rais- 
ing field  to  some  value  //'(each  represented  by  closed 
circle),  then  reducing  field  to  zero.  Dashed  curve 
measured  after  cooling  in  56  kOe  field.  Insert:  reman- 
ent magnetization  vs  //  at  4.2°K. 


In  the  first  set  of  experiments,  the  sample  was 
cooled  in  zero  field  from  above  ~ 100°K  to  various 
temperatures  down  to  4.2°K.  This  procedure  ensured 
that  at  each  temperature  of  measurement  the  sample 
started  in  a magnetically  virgin  state  with  zero  reman- 
ence.  Isothermal  magnetization  curves  were  then  ob- 
tained with  increasing  fields  up  to  56  kOe;  a typical 
selection  of  them  are  shown  in  Fig.  1.  The  initial 
slopes  of  these  curves,  corresponding  to  the  initial 
susceptibilities,  are  plotted  vs  temperatures  in  the  in- 
sert of  the  figure.  The  susceptibilities  clearly  descnbe 
a sharp  maximum  at  ~ 85°K,  in  qualitative  agreement 
with  previous  a.c.  measurements.3  It  should  also  be 
noted  that  at  this  same  temperature  the  magnetization 
curves  as  a whole  undergo  a qualitative  change  in 
shape.  The  curves  above  85°K  are  uniformly  concave- 
downward,  while  the  curves  below  S5°K  start  with  a 
positive  curvature  and  become  concave-downward 
only  at  higher  fields.  With  decreasing  temperature, 
the  inflection  point  of  the  latter  curves  moves  steadily 
upward  in  field,  reaching  ~ 25  kOe  at  4.2°K. 

In  the  second  set  of  experiments,  we  investi- 
gated the  possibility  that  the  anomalous  magnetic  be- 
havior below  85°K  may  be  irreversible.  The  alloy 
sample  was  cooled  in  zero  field  to  4.2°K,  its  field 
raised  to  some  value  H\  and  its  magnetization  meas- 
ured as  the  field  was  reduced  from  H to  zero.  This 


entire  procedure  was  repeated  for  different  values  of 
H\  and  our  results  are  illustrated  by  the  solid  curves 
in  Fig;  2.  Hysteresis  obviously  does  occur  and  gives 
rise  to  an  isothermal  remanence  oR  However,  this  ef- 
fect remains  extremely  small  until  //' is  raised  out  of 
the  initial  linear  region  of  the  virgin  magnetization 
curve.  The  plot  of  oR  vs  //  in  the  insert  of  Fig.  2 
shows  that  the  behavior  is  reversible  below  a threshold 
//'value  of  ~ 15  kOe.  At  higher  values  of  //',  oR  rises 
very  rapidly,  although  its  rate  of  rise  is  decreasing  at 
our  largest  //'value  of  56k0e,  which  suggest  an  ap- 
proach to  saturation.  The  dashed  curve  in  Fig.  2. 
measured  after  cooling  to  4.2°K  in  56  kOe.  gives  3 
remanence  value  of  nearly  6emu/g.  which  is  probably 
also  the  saturation  level  of  the  isothermal  remanence 
at  4.2°K.  The  measurements  of  isothermal  remanence 
were  also  carried  out  at  higher  temperatures,  and  it 
was  found  that  its  threshold  behavior  and  its  approach 
to  saturation  occur  at  //'values  that  decrease  to  zero 
- and  the  remanences  themselves  decrease  to  zero  - 
as  the  temperature  approaches  85°K..  Above  this  tem- 
perature, the  magnetization  is  reversible  at  all  fields. 

Thus,  the  temperature  ( Tm  = 85°K)  of  the  low- 
field  susceptibility  peak  marks  a transition  in  the 
nature  of  the  isothermal  magnetization  processes. 
When  cooled  to  below  Tm  in  zero  field,  the  alloy  ap- 
pears to  sink  into  a low -magnetization  state,  from 
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wlucli  it  can  be  lifted  isotliertnally  but  irreversibly  by 
the  application  ot  a field  above  a certain  threshold 
value.  All  aspects  of  this  behavior  are  consistent  with 
the  onset  of  local  anisotropy  at  Tm  and  its  gradual 
strengthening  at  lower  temperatures.  However,  anis- 
otropy cannot  be  solely  responsible  for  the  low- 
niagncti/ation  state  if  substantial  exchange  forces  are 
present,  which  is  bound  to  be  the  case  especially  in 
this  concentrated  alloy.  If  the  mictomagnetic  picture 
is  valid,  exchange  interactions  in  the  alloy  produce  an 
inhomogeneous  magnetic  state  in  which  all  the  Mn 
moments  tend  to  cancel  out  vectorially.  An  interest- 
ing question  arises  at  this  juncture:  does  the  micto- 
magnetic state  form  simultaneously  witli  the  onset  of 


anisotropy  at  Tm  or  is  it  produced  at  some  higher 
temperature  by  isotropic  exchange  and  is  then  simply 
frozen  in  by  the  anisotropy  at  and  below  Tm1  The 
latter  possibility  appears  to  be  supported  by  recent 
neutron  scattering  experiments  on  a Cu  -Mn  alloy  of 
siinnar  composition,  which  indicate  that  the  magnetic 
correlations  in  the  alloy  arc  essentially  unchanged 
from  4.2  up  to  78°K.5  Further  experiments  are  re- 
quited. however,  before  this  question  can  be  answered 
(or  even  asked)  unambiguously.  Moreover,  the  same 
question  probably  also  applies  to  other  spin-glass 
alloy  systems  such  as  Au-Fe,  whose  isothermal  rem- 
anence  is  found  to  disappears!  the  temperature  of 
the  low-field  susceptibility  peak.6 
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Biquadratic  exchange  coupling  and  the  magnetic  properties  of  PrAg  t 

T.  O.  Bran,  J.  S Kouvel,*  and  G.  H.  Lander 
Arfonnt  National  Laboratory.  Argonnt.  Illinois  60439 
(Received  !5  January  1976) 

Bulk  magnetic  and  neutron  difTraction  studies  show  that  the  cubic  (CsCl-type)  compound  PrAg  becomes 
antiferromagnelic  at  T,  a:  Il’K,  below  which  it  undergoes  a spin-flop  transition  at  a critical  Held  II,  of  - 5 
kOe.  However,  the  temperature  dependence  of  the  initial  susceptibility  below  L„  and  the  field  dependence  of 
the  magnetization  .Vf  above  H.  are  both  anomalous.  The  measured  .Vf  vs  II  curves  for  PrAg  above  7\  and 
those  for  the  related  pseudobinary  Pra)Lx,,Ag  (which  remains  paramagnetic  down  to  4 2"K)  are  tested 
against  the  isotherms  of  ,Vf  vs  effective  field  calculated  for  the  crystal-field  states  of  Pr3*  in  these  compounds 
This  analysis  reveals  that  the  net  exchange  field  H„ch  on  a Pr  atom  in  either  compound  does  not  vary  linearly 
with  .Vf  but  also  contains  a substantial  negative  term  in  Vf  1 It  is  shown  that  this  higher-order  //,«„( .Vf) 
dependence  can  be  generated  by  a positive  biquadratic  exchange  term  in  the  spin  Hamiltonian  of  the 
paramagnetic  system.  Moreover,  When  the  same  biquadratic  exchange  term  is  included  in  a two-sublattice 
model  for  the  ordered  magnetic  state,  a consistent  explanation  is  obtained  for  the  anomalous  properties  of 
PrAg  below  7V  Possible  mechanisms  for  the  origin  of  biquadratic  exchange  in  these  materials  are  discussed. 


I.  INTRODUCTION 

The  cubic  (CsCl-type)  intermetallic  compound 
PrAg,  according  to  recent  neutron-diffraction 
work,1  orders  antiferromagnetically  at  -14°K.  In 
its  collinear  two-sublattice  structure,  the  ordered 
Pr  moments  (~2.1nfl)  in  adjacent  ferromagnetic 
(110)  planes  are  antiparallel  and  lie  within  (001) 
planes.  This  antiferromagnetic  structure  is  identi- 
cal, except  for  the  orientation  of  the  moment  axis, 
to  the  (trtrO)  structures  deduced  from  similar 
studies  of  NdAg,J  TbAg,3  DvAg,8  HoAg,5  and  ErAg,8 
the  latter  two  also  exhibiting  modulated  variations. 
However,  in  two  separate  reports7,0  on  the  bulk 
magnetic  properties  of  this  family  of  isomorphous 
rare-earth-silver  compounds,  it  is  claimed  that 
in  the  case  of  PrAg  the  field  and  temperature  de- 
pendences of  the  magnetization  are  indicative  of 
ferromagnetic  rather  than  antiferromagnetic  or- 
dering. 

This  apparent  inconsistency  prompted  us  to 
undertake  a more  detailed  investigation  of  the 
bulk  magnetic  properties  of  PrAg.  At  low  fields 
C H),  our  magnetization  (A/)  versus  temperature 
data  for  a polvcrystalline  sample  show  a peak  at 
~ 1 1 ’K.  clearly  indicating  the  onset  of  antiferro- 
magnetism. although  there  is  also  an  anomalous 
upturn  in  A/  at  lower  temperatures.  Furthermore, 
the  measured  A/ vs  H curves  below  11°K  exhibit 
an  inflection  point  at  a critical  field  (//,.),  which 
suggests  a spin- flop  transition  to  a canted  state. 

But  again  the  total  behavior  is  anomalous,  in  that 
the  M vs  H curves  at  fields  well  above  hc  are  not 
linear  but  distinctly  concave  downward.  This  high- 
field  behavior,  by  itself,  would  appear  to  agree 
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with  the  previous  reports7’8  that  PrAg  is  ferro- 
magnetic. However,  from  subsequent  neutron- 
diffraction  experiments,  we  have  found  that  PrAg 
does  indeed  undergo  a field-induced  transition  to 
a canted  state.  The  nonlinear  A/  vs  If  behavior  of 
the  spin- flopped  state  thus  remains  as  a real 
anomaly.  Our  experimental  results  for  PrAg, 
which  were  recently  reported  in  brief,9  as  well  as 
those  for  the  related  pseudobinary  compound 
Pr3  5La0-,Ag  (where  the  La  substitution  suppresses 
any  magnetic  ordering  down  to  4.2°K,  at  least), 
are  fully  described  in  Sec.  II  of  this  paper. 

Theoretical  analyses  of  our  results  are  presented 
in  Sec.  III.  Starting  with  the  crystal-field  states 
for  Pr3*  in  PrAg  (where  the  ground  state  is  a mag- 
netic T,  triplet),  as  determined  recently  from  in- 
elastic-neutron-scattering measurements,10  we 
calculate  various  isotherms  of  A/  vs  II c(f , the  total 
effective  field.  Comparing  our  low- temperature 
paramagnetic  A/  vs  H data  against  the  correspond- 
ing A/  vs  Ht„  isotherms,  we  find  that  in  both  PrAg 
and  Pr0  5La0  5Ag  the  net  exchange  field  {fl,^„  = Utt, 

- H)  does  not  depend  simply  linearly  on  A/,  but 
contains  a substantial  negative  term  in  A/3.  We 
proceed  to  show  that  a higher-order  //c<ch(.U)  de- 
pendence of  this  kind  can  be  generated  by  a bi- 
quadratic exchange  term  in  the  effective  spin  Ham- 
iltonian of  the  paramagnetic  system.  Further- 
more. by  the  inclusion  of  the  same  biquadratic  ex- 
change term  in  a two-sublattice  model  for  the 
antiferromagnetic  state,  a consistent  semiquantita- 
tive  explanation  is  obtained  for  the  various  anom- 
alous magnetic  properties  of  PrAg  below  its  Neel 
point. 

In  Sec.  IV  several  possible  mechanisms  for  the 
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origin  of  biquadratic  exchange  in  PrAg  and 
Pr^La^Ag  are  discussed  with  reference  to  pre- 
vious related  work  on  rare-earth  materials. 

II.  EXPERIMENTAL  RESULTS 
A.  Magnetization  and  susceptibility 

For  our  magnetization  measurements,  polycrys- 
talline samples  of  PrAg  and  Pr^jLa^jAg  in  the 
form  of  thin  platelets  weighing  0.13  and  0.11  g,  re- 
spectively, were  cut  from  arc- melted  buttons  and 
annealed  at  700 °C  for  4 h under  argon.  The  same 
annealing  treatment  previously  given  to  larger 
samples  taken  from  the  same  buttons  for  neutron- 
diffraction  study  was  found  to  produce  a single 
cubic  phase  of  appropriate  lattice  spacing.1  The 
sample  preparation  and  their  subsequent  storage 
and  measurement  were  all  carried  out  in  an  inert 
atmosphere  to  avoid  excessive  oxidation.  The 
sample  magnetizations  were  measured  in  a vibrat- 
ing-sample  magnetometer  in  fields  up  to  56  kOe 
at  temperatures  down  to  4.2<1K. 

The  magnetization  (A/,  in  na  per  Pr  atom)  of 
PrAg  is  plotted  in  Fig.  1(a)  against  temperature 
for  various  values  of  the  applied  field  (//);  demag- 
netization corrections  to  //  were  negligible.  In- 
cluded is  the  remanent  magnetization  for  H re- 
duced from  56  kOe  to  zero,  which  is  measurable 
at  low  temperatures  but  extremely  small 
(?0.02us/Pr)  relative  to  the  magnetizations  at- 
tained in  moderate  fields.  A more  prominent  fea- 
ture of  the  curves  for  low  fields  (1  and  3.5  kOe)  is 
the  peak  at  ~11°K.  which  is  typical  of  a polycrys- 
talline antiferromagnet  at  its  Neel  point  (7"v). 
However,  in  that  regard,  the  rise  of  these  curves 
at  lower  temperatures,  corresponding  to  a rise 
in  the  initial  susceptibility,  is  quite  anomalous. 

At  higher  fields,  the  peak  in  the  curves  is  replaced 
by  an  abrupt  change  of  slope,  which  occurs  at  a 
slowly  decreasing  temperature  with  increasing 
field;  below  this  temperature,  the  magnetization 
at  high  fields  is  almost  constant.  This  high-field 
behavior  suggests  that  the  material  has  undergone 
a spin-flop  transition  to  a canted  state. 

Some  of  the  same  data  on  PrAg  are  plotted  in 
Fig.  1(b)  as  isotherms  of  M vs  H . The  curves  for 
12  and  14°K,  in  their  concave-downward  curvature 
at  all  fields,  are  characteristic  of  all  our  data 
above  T„.  Contrastingly,  the  curves  for  4.2  and 
10°K  both  exhibit  an  inflection  point  at  -4.5  kOe, 
which  may  be  taken  as  the  average  critical  field 
Ulc)  for  the  spin-flop  transition  in  this  polycrys- 
talline sample.  However,  above  If c , the  magne- 
tization at  4.2°K  does  not  show  the  linear  increase 
with  field  that  is  normally  expected  of  a spin- 
flopped  state  at  low  temperatures.  Instead,  Af  ap- 
pears to  vary  as  the  cube  root  of  II.  This  unusual 


FIG.  1.  Magnetization  of  PrAg,  in  per  Pr  atom, 
as  (a)  a function  of  temperature  at  various  fields  and 
(b)  a function  of  field  at  various  temperatures. 
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FIG.  2.  .VI2  vs  U/M  for  PrAg  at  4.2 °K,  where  M is 
the  magnetization  (vB  per  Pr  atom)  and  H the  magnetic 
field  (kOe). 


cube-root  dependence  is  brought  out  more  clearly 
in  Fig.  2,  where  the  same  data  are  plotted  as  M2 
vs  ll/M.  Above  a threshold  field,  the  data  in  this 
plot  follow  close  to  a straight  line  that  extrapolates 
back  to  the  origin.  Nevertheless,  despite  its 
anomalous  appearance,  this  high- field  behavior  is 
associated  with  a canted  magnetic  state,  as  our 
neutron-diffraction  study  of  PrAg  in  high  fields 
(described  below)  has  demonstrated. 

Our  low-temperature  magnetization  data  on 
Pr^jLa^Ag  are  plotted  as  isotherms  of  M vs  // 
in  Fig.  3.  From  the  linearity  of  these  curves  at 
low  fields  and  the  smooth  monotonic  changes  with 
field  and  temperature,  it  is  evident  that  this  mate- 
rial remains  paramagnetic  down  to  4.2°  K.  It  has 
therefore  allowed  us  the  advantage  of  extending 
our  analysis  of  the  paramagnetic  state  (Sec.  MA) 
to  temperatures  well  below  the  T „ of  PrAg  and. 
thus,  to  higher  levels  of  induced  moment  per  Pr 
atom. 

In  Fig.  4,  our  low-field  results  for  PrAg  and 
Pr0  sLa^Ag  are  exhibited  as  plots  of  the  inverse 
initial  susceptibility  (x0*')  versus  temperature  ( T). 
The  two  curves  are  essentially  identical  in  shape 
with  a constant  vertical  separation  that  presumably 
arises  from  a difference  in  the  net  exchange  cou- 


FIG.  3.  Magnetization  of  Pr5.5La5.5Ag,  in  p s per  Pr 
atom,  as  a function  of  field  at  various  temperatures. 


FIG.  4.  Inverse  initial  susceptibility  versus  tempera- 
ture for  PrAg  and  Pr0_5La0-SAg.  The  dashed  curve  was 
calculated  for  the  crystal -field  parameters,  IV  =-0.43 
meV,  t =-0.8;  the  dotted  curve  represents  the  r5  ground 
state  alone. 
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pling.  They  both  decrease  in  slope  with  increasim' 
T up  to  ~25°K,  above  which  they  are  linear  and 
obey  the  Curie-Weiss  law,  xi T’ = (T  =0)/Ccw.  with 
C’cw  =0.265  K pflAOe,  corresponding  to  an  effective 
paramagnetic  moment  per  Pr  atom  (pA)  of  3.44(ia. 
The  values  of  0 are  6°K  for  PrAg  and  0°K  for 
Pr0  5 La0>s Ag.  Our  results  for  PrAg  agree  quite 
closely  with  those  of  Walline  and  Wallace7  (p, 

= 3.44(1  a,  O = 2°K).  but  less  closely  with  those  of 
Pierre  and  Pauthenet8  (/i>  = 2.97(ib,  0 = 14°K). 

The  values  of  should  be  compared  with  the  theo- 
retical value  Of  3.58(is  for  the  free  Pr7*  ion. 
which  is  represented  in  Fig.  4 by  the  high-temper- 
ature slope  of  the  dashed  curve  calculated  for  the 
crystal- field  parameters  to  be  discussed  later. 

The  discrepancy  with  the  high-temperature  slopes 
of  our  experimental  curves  conceivably  arises 
from  changes  in  the  exchange  coupling  produced 
by  thermal  lattice  expansion. 

B.  Neutron  diffraction 

To  investigate  specifically  the  possibility  that 
PrAg  undergoes  a field  induced  spin-flop  tran- 
sition. we  carried  out  an  elastic-neutron-diffrac- 
tion study  of  this  material  at  the  CP- 5 reactor  of 
the  AML.  The  sample  was  a 2-g  compact  of  pow- 
der that  had  been  used  for  previous  neutron-dif- 
fraction work.1  All  the  measurements  were  made 
with  the  sample  at  ~ 5 Tv  and  in  a magnetic  field 
of  up  to  60  kOe  applied  perpendicular  to  the  scat- 
tering vector  (k).  Thus,  for  each  ( hkl ) reflection 
investigated,  the  only  crystallites  of  the  sample 
that  contributed  to  the  measured  intensity  were 
those  whose  (hkl)  planes  contained  the  direction  of 
the  applied  field.  The  reflections  studied  were  the 
superlattice  (jjO)  and  (jjl)  and  the  fundamental 
(100). 

The  geometry  for  the  (j^O)  reflection  measure- 
ment is  illustrated  in  Fig.  5.  A scattering  crystal- 
lite is  represented  by  its  chemical  unit  cell  and  the 
Pr  moments  are  shown  in  the  (.ttO)  configuration, 
lying  within  (001)  planes,  as  determined  previous- 
ly in  zero  field.1  Moreover,  although  the  orien- 
tation of  the  spin  axis  within  the  (001)  planes  is 
unknown,  we  assume  for  this  illustration  that  the 
magnetocrystalline  easy  axes  are  the  [110]  and 
[110) , corresponding  respectively  to  cases  A and 
R in  the  figure.  These  two  cases  will  thus  rep- 
resent the  two  types  of  antiferromagnetic  domains 
in  the  crystallite.  In  this  zero-field  situation,  the 
B domains  with  their  spin  axes  normal  to  x con- 
tribute fully  to  the  measured  (}j0)  reflection, 
whereas  the  A domains  with  their  spin  axes  paral- 
lel to  x make  zero  contribution.  When  the  ex- 
ternal field  (H)  is  applied,  it  lies  in  some  direc- 
tion within  the  (110)  plane  of  the  crystallite,  as 


indicated  in  the  figure.  The  field  therefore  has 
a component  parallel  to  the  spin  axes  of  the  11  do- 
mains, which  favors  a reorientation  (or  flop)  of 
the  moments  into  the  plane  perpendicular  to  H. 
Under  these  circumstances,  the  spin  flop  would 
cause  some  reduction  of  the  observed  intensity. 

If  there  is,  in  addition,  a strong  anisotropy  that 
holds  the  moments  in  the  (001)  plane,  the  spin 
flco  induced  by  II  would  effectively  transform  the 
IS  domains  into  A domains,  and  the  observed  in- 
tensity would  reduce  to  zero.  Since  the  component 
of  H parallel  to  the  spin  axes  of  the  IS  domains 
will  differ  for  different  crystallites,  the  latter  pro- 
cess can  be  expected  to  proceed  rapidly  and  then 
more  gradually  throughout  the  sample  as  the  field 
is  increased  beyond  a critical  value.  A similar 
reduction  can  be  predicted  for  the  (((0)  intensity 
with  increasing  field  if  the  easy  axes  are  taken  to 
be  the  [100]  and  [010]  and,  again,  if  the  anisotropy 
out  of  the  (001)  plane  is  very  strong. 

The  intensities  of  the  three  measured  reflections 
at  various  fields  are  shown  in  Fig.  6.  For  the  (^0) 
reflection  with  increasing  field,  the  intensity 
drops  rapidly  between  ~ 4 and  6 kOe  and  then  dimin- 
ishes more  slowly  until  it  has  vanished  at  ~ 60kOe. 
Qualitatively,  this  behavior  is  entirely  consistent 
with  that  described  above  for  the  occurence  of  a 
spin-flop  transition.  Moreover,  the  critical  field 
for  this  transition  appears  to  be  - 5 kOe,  in  good 
agreement  with  the  Hc  value  derived  from  our  mag- 
netization data.  For  the  (£jl)  reflection,  the  mea- 
sured intensity  shows  little  change  until  it  de- 
creases slowly  with  increasing  field  above  - 20  kOe. 
In  this  case,  if  the  spin  axis  of  each  crystallite 
stays  within  the  (001)  plane  following  the  spin-flop 
transition,  as  described  above,  there  should  be 
almost  no  change  of  intensity.  The  intensity  de- 
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FIG.  5.  Schematic  representation  of  neutron-diffrac- 
tion measurements  of  the  < *. i 0)  reflection  from  a PrAg 
crystallite  in  an  external  field  in  the  (110)  plane.  The 
solid  and  dashed  spin  vectors  in  the  (001)  plane  cor- 
respond to  two  types  of  antiferromagnetic  domains. 
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crease  seen  at  higher  fields  can  be  readily  as- 
cribed to  an  increased  canting  of  the  moments  and 
to  Uie  consequent  reduction  of  the  antiferromag- 
netic spin  component,  the  (|^1)  intensity  depend- 
ing on  the  square  of  this  component.  Hence,  the 
fact  that  the  \ 1)  intensity  is  still  appreciable  at 
60  kOe  indicates  that  the  spin  canting  persists  up 
to  much  higher  fields  before  ferromagnetic  align- 
ment is  ultimately  achieved.  The  intensity  of  the 
(100)  reflection,  which  rises  monotonically  with 
increasing  field,  contains  a small  field-independ- 
ent nuclear  part  plus  a magnetic  part  that  varies 
as  the  square  of  the  ferromagnetic  spin  component. 
In  Fig.  6,  we  have  allowed  for  the  nuclear  inten- 
sity and  fitted  the  data  points  with  a curve  derived 
quantitatively  from  the  .Vf  vs  H curve  for  4.2  TC  in 
Fig.  1(b).  The  fit  is  clearly  excellent.  In  essence, 
this  latter  result  confirms  the  cube-root M vs  U 
behavior  of  PrAg  shown  in  Fig.  2.  It  thus  estab- 
lishes, when  combined  with  our  results  for  the 
(||0)  and  (|jl)  reflections,  that  this  unusual  be- 
havior is  exhibited  when  the  material  is  in  a canted 
magnetic  state  following  a spin-flop  transition. 


FIG.  6.  Neutron  intensities  of  various  reflections  from 
PrAg  at  ~5’K,  as  functions  of  the  applied  field  normal 
to  the  scattering  vector.  The  (100)  curve  was  derived 
from  bulk  magnetization  data. 


III.  THEORETICAL  ANALYSES 


A.  I'amn.-ignclic  stale 

Inelastic-neutron-scattering  work10  has  recently 
established  the  crystal-field  level  scheme  for  the 
Pr*  ion  in  both  PrAg  and  Pr0,.  La0_,  Ag.  As  shown 
in  Fig.  7,  the  ground-state  multiplct  of  the  free 
ion  =4,  g =|)  is  split  by  the  cubic  crystal 
field  into  four  levels  which,  with  increasing  energy, 
are  a r5  triplet,  a r3  doublet,  a r.,  triplet,  and 
a r3  singlet.  In  the  notation  of  Lea  et  al.,“  the 
sequence  and  separation  of  these  levels  is  given 
by  iF  = -0.43  meV,  x=-0.8.  This  sequence  is 
the  reverse  of  that  for  the  cubic  (NaCl-type)  Pr 
chalcogenides  and  pnictides,  which  are  singlet- 
ground-state  systems.12'13 

The  ground  state  of  PrAg  and  Pr^,  La0 ,,  Ag  being 
a magnetic  triplet,  it  follows  that  a Pr  atom  (act- 
ing as  a tripositive  ion)  in  these  compounds  will 
be  increasingly  polarizable  by  an  effective  field 
(//„,)  as  the  temperature  decreases  to  zero.  In 
fact,  at  T = 0,  the  initial  effective  susceptibility 
(XorPo  = (M///Cfr)0  diverge.  This  is  shown  in 
Fig.  4,  where  the  calculated  dashed  curve  of 
(x.fr)o 1 vs  T emerges  from  the  origin.  Moreover, 
the  slope  of  this  curve  at  T = 0 is  determined 
solely  by  the  r5  ground-state  triplet  and  is  rep- 
resented by  the  dotted  line  in  Fig.  4.  Its  value 
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corresponds  to  a paramagnetic  moment  (p,)  of 
2.82/ia,  since  the  T,  triplet  alone  constitutes  a 
system  of  pseudospin  6'  = 1,  g =2.  Note  that  the 
calculated  dashed  curve  in  Fig.  4,  in  its  decrease 
of  Slope  with  increasing  temperature,  matches 
the  shape  of  the  experimental  curves  up  to  ~ 10  °K; 
the  discrepancy  at  higher  temperacure  was  com- 
mented on  earlier. 

As  standards  of  comparison  for  the  entire  para- 
magnetic .1/  vs  11  behavior  of  the  two  compounds 
(not  merely  their  initial  susceptibilities),  we 
calculated  various  isotherms  of  M vs  Uc„  on  tne 
basis  of  the  crystal-field  levels  shown  in  Fig,  7. 
The  effective  field  Ht,<  was  taken  to  be  parallel, 
in  turn,  to  each  of  the  principal  crystallographic 
axes.  The  results  of  these  calculations  are  rep- 
resented conveniently  as  isotherms  of  .1/  vs  /£,,  /T 
in  Fig.  8.  The  curve  labeled  0 °K  corresponds  to 
the  isotropic  response  of  the  F,  ground-state  trip- 
let alone;  it  is  a Brillouin  function  and  saturates 
at  2,ifl.  The  curve  labeled  *>  gives  the  isotropic 
high-temperature  response  of  all  nine  states;  its 
saturation  level  is  3.2^s.  These  two  extreme 
curves  bracket  the  curves  for  all  intermediate 
(normal)  temperatures  for  all  directions  and 


Heff/T  (kOe/K) 

FIG.  8.  Magnetization  (in  Ug  per  Pr  atom)  vs  W,ir/T 
calculated  for  PrAg  and  Pr0_jDi0  sAg  at  various  tempera- 
tures (7",  in  °K)  for  effective  fields  (//,„,  in  koe)  paral- 
lel to  (100)  or  (100) , 


n 

strengths  of  II,,,.  At  each  temperature.  Iho  initial 
slopes  of  the  curves  are  the  same  for  all  direc- 
tions of  He,  i.e.,  the  initial  susceptibilities  are 
isotropic,  as  required  by  the  cubic  symmetry  of 
the  structure.  However,  with  increasing  field, 
the  curves  for  H,,f  parallel  to  (111)  diverge  from 
and  lie  consistently  above  the  corresponding  iso- 
therms for  H,(r  along  ( 100) ; the  curves  for  Hc„ 
along  (110)  are  intermediate  and  are  not  shown  in 
Fig.  8.  At  the  highest  magnetizations  (-  1.5/ia/Pr) 
reached  by  our  paramagnetic  data  for  the  two  com- 
pounds, the  anisotropy  of  these  curves  is  fairly 
small;  for  our  data  analysis,  we  specifically  used 
the  curves  for  H„,  along  (100) . 

In  the  analysis  of  our  paramagnetic  data  for 
PrAg  and  Pr^La^Ag,  the  net  effect  of  all  the 
exchange  interactions  on  each  Pr  atom  is  consider- 
ed to  be  representable  by  an  exchange  field 
which  is  collinear  with  the  externally  applied  field 
H.  With  the  direction  of  H as  reference,  the  mag- 
nitude of  the  total  effective  field  is  given  by  tlcii 
= //„«„+//,  where  // , ,,  may  be  positive  or  negative. 
Hence,  within  this  effective-field  approximation, 
we  can  consider  an  experimental Mill,  T)  data 
point,  refer  to  the  appropriate  isotherm  among  the 
calculated  curves  represented  in  Fig.  8 and  find 
the  value  of  corresponding  to  the  measured  At. 
and  tnus  determine  the  exchange  field.  //c!( 

-11 , for  the  particular  M and  T.  This  procedure 
was  followed  with  all  our  low-temperature  paramag- 
netic data  on  the  two  compounds,  including  the 
data  on  PrAg  for  11  -s  7 kOe  at  10  °K,  which  from 
Fig.  1(a)  appear  to  be  just  within  the  paramag- 
netic regime. 

The  results  of  this  data  analysis  are  presented 
in  Fig.  9(a)  as  isothermal  plots  of  M 2 vs  ffeVvtl/.Vf. 

If  simple  mean-field  theory  is  valid  and  is 
proportional  to  .1/  (and  not  otherwise  dependent 
on  T),  the  reduced  data  points  in  this  form  of 
representation  can  be  expected  to  follow  a single 
vertiele  straight  line  for  each  material.  Clearly, 
the  results  in  Fig.  9(a)  do  not  bear  out  this  ex- 
pectation. Although  the  isotherms  for  each  com- 
pound do  cluster  fairly  closely  together  and  they 
do  form  reasonably  straight  lines,  the  lines  have 
a pronounced  tilt  corresponding  to  decreases  of 
f/„lh/Af  with  increasing  A/2.  At  high  U2,  the  iso- 
therms appear  to  be  deviating  from  linearity  to- 
wards even  smaller  values  of  ffcub  If.  In  fact, 
in  the  case  of  Pr„,sLa, . ,Ag  at  4.2  K, 
ultimately  changes  sign  from  positive  to  negative — 
which  underscores  the  enormous  magnitude  of 
this  phenomenon.  These  results  can  be  expressed 
in  summary  as  follows: 

(l) 

or.  equivalently, 
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FIG.  9.  Isotherms  of  .V/2  vs  (H.ff  -H)/M  derived  from 
data  on  PrAg  and  Prn  sLa5-5Ag  on  the  basis  of  (a)  the 
crystal -field  parameters,  If  = — 0.43  meV,  x = — 0.3, 
and  (b)  the  Fs  ground  state  alone. 

, (1') 

where  a is  the  usual  exchange  coefficient  and  a' 
is  the  coefficient  associated  with  some  higher- 
order  process.  Terms  of  still  higher  order  in 
M are  implied  in  these  expressions  in  order  to 


allow  for  the  nonlinearity  of  the  isotherms  in 
Fig.  9(a)  at  high  A/2;  they  should  be  such  that 
odd  symmetry  between  //xch  and  M is  preserved. 
The  values  of  A and  A'  derived  from  the  data  in 
Fig  9(a)  are  listed  in  Table  I.  For  each  com- 
pound, the  two  coefficients  show  little  variation 
with  temperature. 

For  comparison,  we  carried  out  an  analogous 
treatment  of  the  same  data,  using  exclusively 
the  .U(//c„,  T ) characteristics  of  the  F5  ground - 
state  triplet  (i.e.,  the  0 ’K  curve  in  Fig.  8).  The 
results  for  ,U2  vs  //eNt„/.V/  are  plotted  in  Fig.  9(b). 
The  isotherms  again  form  lines  of  decidedly  neg- 
ative slope.  However,  their  values  at  zero  M 
and  their  slopes,  corresponding,  respectively, 
to  A and  (A')"1  in  Eq.  (1),  vary  much  more  with 
temperature  than  those  of  their  counterparts  in 
Fig.  9(a).  The  A and  A'  values  are  listed  in  Table 
I and  bear  out  this  comparison.  Thus,  if  the 
constancy  cf  A and  A'  with  temperature  is  a mean- 
ingful criterion,  it  appears  distinctly  preferable 
in  this  type  of  analysis  to  use  the  M(Ht„ , T)  curves 
computed  for  the  full  spectrum  of  crystal-field 
levels  rather  than  the  curves  for  the  T,  triplet 
alone.  Nevertheless,  since  the  results  in  Fig. 

9(al  are  quite  similar  in  basic  character  to  those 
in  Fig.  9(b),  it  is  not  unreasonable  at  these  low 
temperatures  to  consider  only  the  ground-state 
triplet  and,  in  doing  so,  simplify  enormously 
any  theoretical  investigation  of  the  magnetic  pro- 
perties of  the  two  compounds. 

The  higher-order  dependence  of  H ’.yil  on  .U.  as 
revealed  by  our  data  analysis,  presumably  arises 
from  an  equivalently  higher-order  component  of 
the  effective  exchange  coupling.  Specifically,  the 
,U3  term  in  Eq.  (1')  suggests  that  the  spin  Ham- 
iltonian for  these  materials  contains  an  exchange 
term  that  is  quartic  in  the  spin  variables.  The 
simplest  of  such  terms  has  the  biquadratic  form 
(S,  - S,  )2  for  two  interacting  spins.  In  keeping  with 
the  fact  that  all  our  experimental  data  are  for 
polycrystalline  samples,  we  will  restrict  our 


TABLE  I.  Bilinear  and  biquadratic  exchange  coefficients. 

PrAg 


r (°k) 

If  = -0.43 
A (kOe/i i8) 

meV,  t = -0.3 

X (kC  )o/ IX  g) 

r5  is=i 

A (kOe/  us) 

X'  (kOe/^g) 

10 

36.  G 

—0.3 

45.  0 

-4.0 

12 

33.5 

— 0.4 

50.3 

-3.7 

14 

39. 1 

—0.4 

53.0 

-3.2 

Pr0s 

I^VsAg 

4 

12.4 

-5.3 

14.4 

-5.1 

0 

16.6 

-5.3 

20.1 

-4.5 

10 

19.2 

— 5.3 

27.0 

-3.5 
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considerations  to  this  isotropic  biquadratic  term 
plus  an  isotropic  bilinear  exchange  term  of  the 
usual  form  S,  • S, . Moreover,  in  the  approxima- 
tion that  PrAg  and  Pr^,  5La„.  sAg  are  characterized 
adequately  by  their  crystal-field  ground-state 
triplet,  they  will  be  regarded  as  simple  isotropic 
systems  of  effective  spin  3 = 1 and  xr  - 2.  Extending 
from  this  particular  case  to  any  5,  we  express 
die  quantum-mechanical  Hamiltonian  of  the  itn 
spin  of  the  system  as 

= £ [«/„§.■§, +J;,(S(«S1)»],  (2) 

j «i 

where  H is  the  external  field,  and  Ju  and  J\,  are, 
respectively,  the  coefficients  for  the  bilinear 
and  biquadratic  exchange  interactions  with  all 
other  (jth)  spins. 

The  mean-field  approximation  to  the  spin  Ham- 
iltonian of  Eq.  (2)  can  be  written  as 

= ~ iLt  ( ^ (S  ct  ii \^ia  ) )S|a 

i " i \ a 

aSj  g)5j  aSlB\  , (2  ) 

a.fl  / 

where  a and  3 represent  Cartesian  coordinates 
(x.y,  z)  permuted  in  the  summations,  and  ( ) 
corresponds  to  a thermal  average.  This  mean- 
field  Hamiltonian  has  been  treated  previously  by 
Lines  and  Jones14  for  the  case  of  a collinear  anti- 
ferromagnet  in  zero  H.  Our  problem  here  being 
the  paramagnetic  state  in  nonzero  H,  we  will 
apply  the  same  treatment  but  with  H placed  along 
the  z axis  and  with  the  ] summations  making  no 
distinction  between  magnetic  sublattices.  Thus, 
following  Lines  and  Jones,  we  obtain  from  Eq.  (2') 

~ j — it 

i • I 

h;.  )| . 


the  Stevens  quadrupole  operator,  0°  = 3S;  - S(S  + 1), 
the  last  two  terms  on  the  right  in  Eq.  (4)  become 
-K[<0;)0^-2S(S  + 1)Ss(],  thus  showing  the  well- 
known  equivalence  between  biquadratic  and  quad- 
rupole coupling.  These  terms  reduce  to  a con- 
stant for  S = i;  biquadratic  interactions  have  no 
effect  on  the  thermodynamic  properties  of  a spin- 
i system.  Since  3C(  in  Eq.  (4)  is  diagonal  in  Sl£ , 
the  eigenenergies  are 

- ^[3<Si>-S(S+l)]m- 

+ 5^[<Si)-S(S+l),  (5) 

where  in  is  an  eigenvalue  of  Sit.  For  consistency 
within  the  model,  (Sit)  = (S£)  and  (S^)  = (S^), 
and  the  expressions  for  these  thermal  averages 
are 


(S,)=  t6a) 

*»  "I 

and 

JV**.'*1-,  16b) 

where  the  summations  scan  the  m values:  -S, 
-S+l S-  1.  S. 

Although,  in  principle,  a self-consistent  solu- 
tion of  Eqs.  (5),  (6a).  and  f6b)  would  give  \S£'  as 
a function  of  H , and  T . the  complexity  of  the  pro- 
blem under  general  conditions  would  require  that 
the  solution  be  numerical.  To  obtain  an  analytical 
solution,  we  resort  to  an  expansion  effectively  in 
terms  of  J'a/kT . The  details  of  this  procedure 
are  described  in  Appendix  A.  We  show  that  in 
the  equation  of  state  written  as 

(S,)/S*B,UnBSH,tt/kT),  (7) 

where  Bs  is  the  Brillouin  function,  the  total  effec- 
tive field  H't,  can  be  expressed  as  follows: 


(3) 

By  virtue  of  tbe  fact  that  S;,  *3“,  =5;  and  that 

<S-;,>  = <$5,>  = i[S(S+  1)  - (S},)), 

Eq.  (3)  can  be  written  as 

3C(  = -Qrp*f/,+(</0-^XS,>!sl, 

-i^[3<s-;>-s(s+i)js-, 

+ i^«SJ,>-S(S  + l)]S*,  (4) 

where  J’a^H  l,cl’u,  and  the  ; in  < SJt > 

and  (S:j£)  has  been  dropped  because  all  the  mag- 
netic atoms  are  considered  equivalent.  Lf  we  use 


n 

+ Lww)v;(s,)!+-,  (a) 

n 

where  n = 0,l,2 The  coefficients  .4„  and 

are  numerical  functions  of  S (which  arc  zero  for 
any  n if  S = '.)  and  are  given  in  Appendix  A . Equa- 
tion (8)  can  obviously  be  converted  to  an  expres- 
sion for  the  exchange  field  - H t),  which 

for  J’a,  « AT  is  simply 

» (*  o - J* o*Va)3/  + A0.V  ' M 3 + • • • , 


1 


(9) 


13 


BIQUADRATIC  EXCHANGE  coupling  and  the  magnetic. 


5015 


Where  ^0SJo/sz^a>  and  M=gng(S). 

the  average  moment  per  magnetic  atom. 

Since  A. ' represents  the  biquadratic  part  of  the 
spin  Hamiltonian,  Eq.  (9)  shows  that  biquadratic 
interactions  in  a paramegnetic  system  (of  S— 1) 
cause  the  net  exchange  field  to  contain  cubic 
(and  higher-order)  terms  in  the  induced  magnet- 
ization. This  result  clearly  provides  a plausible 
explanation  for  the  paramagnetic  ff,lch  vs  M char- 
acteristics  of  PrAg  and  Pr^La^Ag,  as  expressed 
m Eq.  (1  ).  m relating  \'  and  A„',  the  coef- 
ficients of  A/3  in  Eqs.  (1')  and  (9),  we  note  in 
Table  I that  the  A'  values  for  both  compounds  are 
negative,  and  since  A0  is  positive,  the  correspond- 
ing values  of  A ' are  also  negative.  A negative 
\0,  in  turn,  corresponds  to  a negative  value  for 
the  j summation  of  Ju  in  Eq.  (2),  which  means 
that  predominantly  the  biquadratic  interactions 
favor  quadrature  rather  than  coUinear  spin  align- 
ment. Hence,  in  these  materials,  the  biquadratic 
and  bilinear  interactions  must  be  in  competition 
since  the  latter  always  favor  a coUinear  spin 
alignment. 

In  the  theoretical  problem  of  the  antiferromag- 
netic properties  of  PrAg  in  an  external  field 
which  is  discussed  in  Sec.  UIB,  the  mean-field 
Hamiltonian  analogous  to  Eq.  (4)  would  be  much 
too  complicated  to  allow  analytical  solution.  How- 
ever. by  following  a simplified  mean-field  ap- 
proach. we  manage  to  obtain  analytical  expressions 
that  predict  a physicaUy  reasonable  behavior  for 
each  property  of  interest.  In  order  that  a valid 
analytical  connection  can  later  be  made  between 
the  properties  of  PrAg  above  and  below  T/Il  we 
wUl  now  follow  the  same  simplified  approach  for 
the  paramagnetic  state.  Hence,  we  set  (S,  S ) 

= <Sf'<S*>  in  Eq.  (2'),  which  for  HX=H  =0  becomes 
simply 

K'  = -lniaHt  + J0GJ)Slt  - J'0{St)2S]„  (10) 

where  J0  and  J'  are  defined  as  before.  Here  Xt 
is  diagonal  in  Su  and  thus  gives  immediately 

E«  = -(ZVBH,  + J0(Sl,))>n-J^S,)2>„*  (11) 

for  the  eigenenergies,  where  in  is  an  eigenvalue 
of  Slt.  For  an  analytical  expression  for  H%it , we 
again  resorted  to  an  expansion  effectively  in  terms 
of  </’/  AT.  As  described  in  Appendix  A , when  H 
is  defined  by  Eq.  (7),  it  again  has  the  form  give'n 
by  Eq.  (8),  except  that  is  omitted  from  the  co- 
efficient of  ( S ,).  Thus,  the  exchange  field  can  be 
written  as 

M each  * A„Af  + AqX^AI  3 + • • • f (12) 

where  A 0 and  a^  are  defined  as  before,  and  A„ 
and  higher- order  coefficients  are  given  in  Appen- 


dix A. 

Clearly,  even  in  this  simplified  mean-field  ap- 
proximation, the  biquadratic  interactions  in  a 
paramagnetic  system  generate  cubic  and  higher- 
order  terms  in  M in  the  expression  for  II 
Hence,  similar  to  Eq.  (9),  Eq.  (12)  provides  a 
simple  theoretical  rationale  for  the  paramegneMc 
behavior  of  PrAg  and  Pr0.5La„.5Ag.  Our  qualita- 
tive comments  about  the  negative  values  of  A' 
for  these  materials  are  still  valid.  Quantitatively 
we  can  identify  Eq.  (12)  with  Eq.  (1'),  such  that 
A0=A  and  AJA'  = a',  and  use  the  A and  A'  values  in 
Table  I for  the  two  compounds  considered  as  S = 1 
systems.  Thus,  also  using  AJ  = £ for  S = 1 , we  ob- 
tain for  PrAg  at  10°  K 

A 0 = 45.6  kOe/ pfl,  Ag  = -9.2  kOe/pJ,  (13a) 
and  for  Pr^La^Ag  at  4.2°K 

A0  = 14.4  kOe/pa,  A^  = -10.2  kOe/p3,.  (i3b) 

It  is  obvious  from  these  values  that  the  net  bi- 
quadratic coupling  is  roughly  the  same  in  the  two 
compounds  and,  in  the  case  of  Pr0  ^a^  sAg  is 
comparable  in  strength  to  the  net  bilinear  coupling. 
The  values  of  these  exchange  parameters  for 
PrAg  will  be  discussed  below  with  reference  to 
our  analysis  of  its  antiferromagnetic  behavior. 


tv  ith  regard  to  its  ordered  antiferromagnetic 
state,  we  will  continue  to  characterize  PrAg  by 
its  crystal-field  ground-state  triplet  and  thus 
consider  it  effectively  as  a spin-1  system.  More- 
over, the  Hamiltonian  in  Eq.  (2')  will  be  assumed 
to  remain  valid  below  the  Neel  temperature,  ex- 
cept that  the  interaction  summations  over  neigh- 
boring atoms  will  now  involve  the  spins  on  two 
different  sublattices.  Our  primary  interest  being 
the  canted  spin-flopped  state  (and  its  anomalous 
magnetization  behavior  described  earlier)  we 
set  the  sublattice  moments  along  the  ±z  directions 
and  allow  them  to  deviate  symmetrically  from 
these  directions  in  response  to  an  external  field 
in  the  x direction.  Unfortunately,  under  these  con- 
ditions, the  solution  of  Eq.  (2')  for  the  eigenener- 
gies is  extremely  difficult.  To  make  this  equation 
tractable,  we  set  (SJaS,3)  = (SJa){Sia>,  thus  ne- 
glecting spin  correlations,  and  obtain  the  following 
expressions  for  the  Hamiltonian  of  the  ith  spin  on 
sublattice  A (the  other  sublattice  being  B): 


X * = “ MtS*  - _ jj>x  _ j'xP2  _ 


where 


(14) 


5 i';u.  J[  s * JX 
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An  interchange  of  A and  li  produces  an  equivalent 
expression  for  30^.  Furthermore,  from  symmetry 
we  have  = <S*)»<S,>  and  <S*  >«_  (Sf,)»<S,). 
Hence,  Eq.  (14)  can  be  written  without  super- 
scripts for  either  sublattice  as  follows: 

3C<  = - U'U  *#,+  («/,  *J1)(S')\Si,-(Jl-J3)(S,)Slt 

- (^-^)<SI)<S,)(SlrSu  + SlJSJ,  (15) 

where  subscripts  1 and  2 continue  to  refer  to  the 
interactions  within  and  between  sublattices,  re- 
spectively. 

Though  Eq.  (15)  is  valid  for  any  S,  our  particular 
concern  here  is  the  case  of  S = 1 (g  = 2),  which 
pertains  most  closely  to  PrAg.  Various  aspects 
of  this  case,  corresponding  to  various  thermo- 
dynamic properties  derived  from  the  eigenvalues 
of  the  Hamiltonian  (15),  are  discussed  in  detail 
in  Appendix  B.  In  particular,  the  dependence  of 
the  average  magnetization  M{=giiB{Sx))  on  the 
external  field  H(=HJ  is  found  to  obey  the  following 
relationship  at  7 = 0: 

- 4\f[l-*-4,\,Vi/U1  - A. 3) J A/3  +•  • • , (16) 

where 

lb)‘. 

This  equation  relates  specifically  to  the  spin- 
flopped  state  and,  consequently,  the  reciprocal 
of  the  coefficient  of  .1/  represents  the  perpendicu- 
lar susceptibility  at  low  fields,  i.e., 

X1  = -1/2(.V! -4  Xlp2fl)  (17) 

at  7 = 0.  Aside  from  zero-point  energy  effects, 
the  corresponding  parallel  susceptibility  (\„)  can 
be  expected  to  be  zero.  In  comparison,  the  initial 
susceptibility  at  the  N£el  point, 

x(7v)  = — l/2\j,  (18a) 

and  the  N£el  temperature  itself, 

= (18b) 

as  also  shown  in  Appendix  B. 

Since  \(TS)  is  positive,  X-  has  to  be  negative, 
which  simply  reflects  the  fact  that  the  net  bilinear 
coupling  between  sublattices  is  antiferromagnetic. 

If  XI  is  also  negative,  it  follows  from  Eqs.  (17) 
and  (18a)  and  the  condition,  \„=0,  that  the  powder 
susceptibility,  i(2\1^x  ).  at  7 = 0 may  actually 
exceed  the  susceptibility  at  7V.  Thus,  the  pecu- 
liar low  temperature  increase  that  we  have  ob- 
served in  the  initial  susceptibility  of  a polycrys- 


I 


talline  sample  of  PrAg  can  be  ascribed  to  a nega- 
tive X',  corresponding  to  a net  biquadratic  cou- 
pling between  sublattices  that  favors  quadrature 
spin  alignment.  The  role  of  Xl  in  PrAg  is  made 
even  clearer  when  the  anomalous  magnetization 
versus  field  behavior  of  the  material  in  its  spin- 
flopped  slate  is  interpreted  in  the  light  of  Eq.  (1C). 
As  shown  in  Fig.  2,  the  observed  behavior  at 
4.2°K  is  represented  above  the  critical  field  by 
the  relation,  //  = c.U3,  where  c = 18.5  kOe/Va. 
Hence,  a comparison  with  Eq.  (16),  though  strict- 
ly valid  only  at  7=0,  gives 

XJ-4Xlp-;  = 0,  (19a) 

and 

-4Xl[l  + 4X^g2 /(X,  - Xj)]=  18.5  kOe/Va.  (19b) 

Equation  (19a)  may  be  taken  to  signify  that  in 
PrAg  at  4.2 there  is  a balance  between  the 
negative  biquadratic  and  negative  bilinear  inter- 
actions between  sublattices,  in  their  respective 
favoring  of  quadrature  and  antiparallel  spin  align- 
ment. In  fact,  if  -4Xlnj  were  larger  man  -X,, 
such  that  the  coefficient  of  M in  Eq.  (16)  became 
negative,  this  equation  would  predict  a spontane- 
ous magnetization  arising  from  a canting  of  the 
sublattice  moments  in  zero  external  field  This 
zero-field  canted  state  may  actually  exist  in  PrAg 
at  very  low  temperature,  as  will  be  discussed  be- 
low. 

For  a quantitative  estimate  of  the  exchange  co- 
efficients for  PrAg,  we  start  with  the  results 
given  in  Eq.  (13a)  for  the  paramagnetic  state  close 
to  the  N4el  point.  Since  each  of  x0  and  x£  repre- 
sents ar.  algebraic  sum  of  the  interactions  within 
and  between  sublattices,  these  results  can  be  ex- 
pressed as 


X,  + X-  =45.6  kOe/aa 

(20a) 

Xj>Xl  = -9.2  kOe/V„. 

(20b) 

For  the  individual  bilinear-exchange  coefficients 
(Xi.Xj),  we  substitute  our  experimental  values. 
x(7v)  = 0.0 (0aa/kOe  and  7V  = 11  K,  into  Eqs.  (18a) 
and  (18b)  and  obtain  from  their  combination 

= 54.3  kOe/|xs,  X- = -7.1  kOe/>ia.  (21) 

These  results  give  X,*X3  = 47.2  kOe/ufl,  in  very 
good  agreement  with  Eq.  (20aJ.  It  follows  from 
(21)  that  the  positive  value  of  X,-*.  X-,  which  cor- 
responds in  Fig.  4 to  the  downward  displacement 
of  the  experimental  vs  7 curve  from  the  cal- 
culated (dotted)  curve  for  the  noninteracting  spin- 
1 system  at  7V,  reflects  the  fact  that  the  net  bi- 
linear coupling  m PrAg  is  ferromagnetic  within 
the  sublattices  and  relatively  weakly  antiferro- 
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magnetic  between  them. 

In  determining  the  biquadratic-cxcliange  coef- 
ficients (A[,  xJ.)  individually,  it  might  appear  that 
Eq.  (20b)  can  be  used  in  conjunction  with  Eq.  (19b) 
under  the  assumption  that  the  latter  equation  con- 
tinues to  represent  quantitatively  the  cubic  part  of 
the  M vs  II  curve  for  PrAg  down  to  O'K.  {Pre- 
sumably, Eq.  (19a)  does  not  hold  at  0JK  since  it 
would  require  a continued  balance  between  the  op- 
posing effects  of  \f  and  vf.]  However,  for  the 
and  .V,  values  given  in  (21),  the  left  side  of  Eq. 
(19b)  can  only  reach  a maximum  level  of  15.4 
kOe/Vs,  corresponding  to  A’ = -7.7  kOe/g^.  If 
a of  even  approximately  this  value  were  valid, 
the  coefficient  of  .1/  in  Eq.  (16)  would  be  sizably 
negative,  which  would  signify  a nonzero  spontane* 
ous  magnetization  at  O'K,  as  was  mentioned  ear- 
lier. Unfortunately,  our  measurements  do  not 
extend  below  4.2°K,  and  the  possibility  that  PrAg 
has  a zero-field  canted  state  at  lower  tempera- 
tures can  only  be  offered  as  a speculation  (al- 
though a previous  unpublished  report15  does  pre- 
sent experimental  evidence  for  a zero-field  mag- 
netization at  ~2 "K).  The  only  firm  quantitative 
statement  we  can  make  at  present  about  the  bi- 
quadratic interactions  in  PrAg  is  given  m Eq. 
(20b).  Qualitatively,  however,  we  have  shown 
that  these  higher-order  interactions  are  suffi- 
ciently strong  for  them  to  play  a very  important 
role  in  the  magnetic  properties  of  this  material. 

IV.  DISCUSSION 

Some  of  the  earliest  interest  in  biquadratic  ex- 
change coupling  and  its  effect  on  magnetic  prop- 
erties was  with  regard  to  certain  compounds  of  - 
manganese. In  fact,  the  paper  by  Lines  and 
Jones, u whose  approach  was  followed  in  writing 
the  mean-field  Hamiltonian  of  Eq.  (3),  was  con- 
cerned specifically  with  a-MnS.  From  an  analy- 
sis of  NMR  data  for  the  sublattice  magnetization 
versus  temperature,  they  concluded  that  nearest- 
neighbor  biquadratic  interactions  in  this  antiferro- 
magnetic  compound  play  a fairly  insignificant  role. 
This  is  essentially  what  they  concluded  earlier  re- 
garding MnO,1? 

More  recently,  most  experimental  and  theoreti- 
cal investigations  concerned  with  biquadratic  in- 
teractions have  involved  various  compounds  of 
the  rare  earths.  Particularly  in  the  case  of  ionic 
insulators,  there  are  many  examples  in  which 
these  higher-order  interactions  have  been  shown 
to  be  comparable  in  strength  to  the  normal  bilin- 
ear coupling.18  However,  in  the  case  of  rare- 
earth  intermetallic  compounds,  which  pertains 
more  directly  to  our  study  of  PrAg,  experimental 
evidence  for  the  existence  of  biquadratic  exchange 
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has  only  been  reported  for  a small  number  of  ma- 
terials. Among  these,  only  the  pseudobinary 
(Tb,  Y)Sb  compounds,  which  are  singlet  ground- 
state  systems,  have  been  studied  magnetically 
in  any  manner  analogous  to  that  of  the  present 
work  on  PrAg.  It  appears,  however,  that  although 
the  measured  M vs  I!  isotherms  for  the  paramag- 
netic (Tb,Y)Sb  compounds19  are  consistent  with  a 
higher-order  dependence  of  the  exchange  field  on 
Al,  similar  to  that  of  Eq.  (1*),  the  temperature 
dependence  of  the  sublattice  magnetization  of  the 
antiferromagnetic  compounds20  gives  no  support- 
ive evidence  that  the  exchange  field  varies  other 
than  linearly  with  M. 

To  date,  the  most  conclusive  indication  of  bi- 
quadratic coupling  in  rare-earth  intermetallic 
compounds  has  emerged  from  the  results  of  ul- 
trasonic vibrational  measurements  on  DySb, 21,20 
and  HoSb.22  The  pronounced  elastic-constant 
softening  observed  in  these  materials  just  above 
their  magnetic  (and  structural  transition  tempera- 
tures was  attributed  in  large  part  to  magnetic 
quadrupole-quadrupole  (effectively  biquadratic) 
interactions,  where  the  quadrupoles  are  coupled 
to  the  lattice  strains.  Presumably,  the  same  type 
of  interactions  are  responsible  for  the  elastic- 
constant  softening  seen  in  TmCd  as  this  CsCI- 
structured  compound,  which  does  not  order  mag- 
netically, approaches  its  structural  transition 
temperature  from  above.22  Furthermore,  in  the 
case  of  DySb,  the  first-order  transition  to  an  an- 
tiferromagnetic state  has  been  shown  to  be  con- 
sistent with  a net  biquadratic  coupling  that  is  pos- 
itive and  thus  favors  collinear  arrangement  of  the 
spins.24  It  should  be  noted,  however,  that  DySb 
also  undergoes  a field-induced  transition  to  a 
state  where  the  sublattice  moments  are  in  quadra- 
ture.25 which  suggests  that  at  high  fields  the  ef- 
fective biquadratic  coupling  may  be  negative 
(which,  as  we  have  shown,  is  the  situation  in  PrAg 
at  all  fields).  Pertinent  to  this  question,  an  ex- 
tensive phenomenological  study  has  been  made  of 
the  different  types  of  ordered  states  that  can  oc- 
cur at  various  temperatures  in  a spin-1  system 
due  to  the  coexistence  of  bilinear  and  biquadratic 
coupling  of  either  sign.29  Since  this  study  was  re- 
stricted to  collinear  magnetic  ordering,  it  does 
not  encompass  the  possibility  of  a canted  spin 
state,  whereas  such  a state  is  expressly  allowed 
(even  at  zero  field)  by  our  simplified  Hamiltonian 
analysis  of  the  magnetic  properties  of  PrAg  below 
its  N6el  point. 

In  theoretical  discussions  of  DySb,21,22,24,27  the 
biquadratic  interactions  have  generally  been  as- 
cribed to  either  or  both  of  two  sources.  One  is 
the  mechanism  whereby  the  magnetic  quadrupoles 
of  the  rare-earth  atoms  are  coupled  via  the  co- 
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licrent  dynamical  distortions  of  (he  local  crystal 
fields,  produced  by  the  phonons  in  the  system. 
Tins  "virtual  phonon  exchange”  process,  which 
has  recently  been  reviewed  in  detail, 28  can  Rive 
rise  to  an  effective  biquadratic  coupling  of  con- 
siderable range.  The  other  biquadratic  coupling 
mechanism,  which  also  is  inherently  long  range 
and  is  specifically  relevant  to  intermetallics,  is 
a higher-order  component  of  the  indirect  exchange 
interactions  via  the  conduction  electrons.  Higher- 
order  indirect  exchange  coupling  has  recently  been 
calculated  to  be  of  significant  strength  in  DyZn,29 
which  resembles  PrAg  in  its  CsCl-type  structure. 
Although  it  has  been  impossible  to  distinguish  be- 
tween these  two  mechanisms  experimentally  on 
the  basis  of  available  data,  we  are  planning  to  ac- 
complish this  for  PrAg  by  means  of  detailed  neu- 
tron scattering  measurements  of  the  magnon  and 
phonon  spectra. 
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APPENDIX  A 

Following  the  notation  of  Lines  and  Jones,11 
w define 

(Ala) 

and 

l'  = W;i3(S2;  -S(S -I’ll.  (Alb) 

so  that  Eq.  (5)  can  be  written  simply  as 

= - -V*h  - Yin1  * £0,  \A2) 


where  J - 1 kT  and  ms-S.-Sil 5-1,5, 

Under  the  condition  that  J(A '.m  «■  l'<u2)«l  for  any 
m,  the  exponential  is  expanded  as 

exp[d(.Y«i  *•  FmJ)|  = 

*'o.. . . . . 

where  = 1,  £,  =dA\  f , = JV  » >,  0tc..  and 

Eqs.  (A3a)  and  (A3b)  become 

~ lATa) 

• / « 

and 

tA4b> 

V / V 

where  a,  = in *,  which  is  nonzero  only  for  If 
even. 

From  (A4b)  and  (Alb),  we  obtain  1‘  as  a series 
expansion  in  even  (but  nonzero)  powers  of  .V, 
which  we  then  substitute  into  (A4a!  and  obtain 
as  a series  expansion  in  odd  powers  of  A. 

The  latter  series  is  inverted  to  give 

A*  = T i M (Ji  ~ J *■  (>'5  - P5)(4i,/5  ■*■  • * • , 

(A5) 

which  is  so  written  tliat  the  •, 's  depend  only  on 
J and  6'  and  the  p's  are  functions  of  J’  as  well 
as  ,1  and  5.  Moreover,  we  obtain 

r,\S<)  + •••  = S). 

where  B~‘  is  the  inverse  Brilloum  function. 

Hence,  (A5)  may  be  expressed  as 

41).  (All) 

From  this  form  of  the  equation  of  state,  it  is 
clear  that  the  terms  between  the  brackets  con- 
stitute gpg//,,,  , where  //  ,,  is  the  total  effective 
field  on  each  paramagnetic  atom.  Thus,  front 
the  definition  of  A in  (Ala), 

S^b^  II  ♦ f'fo  “ u o ) (4 

-P,\b,)5  * ••  • . (AT) 

where  it  is  found  that 


where  E0  = vJ'0[\?a)  - SiS  «- 1 )] . Substituting  (A2) 
into  Eqs.  (Sa)  and  (6b)  and  factoring  out  E„,  we 
obtain 


(5,1  exPi  4(A'»n  * Yin1)].  expf  ,J(A'm ♦ Ynr) 


and 


(A3a) 


exP[  4(A'm  - Yin1)]  j £exp|  d(-V//;  * Yin1)} , 


lA3b) 


P,  = y;  AA4J’0yr0,  p,  - £ . . . ; 

q ~ A „ = 7]  h3/  ,.4 , — j for  n ‘1; 

^0  a ITS  by  4 ^-J/l«  P(4>»1).  2#,  = I;2  4 * 1 ), 

Ml  1 -</  i'(i  + l)|  '(1  - if/). 

Note  that  which  is  contained  in  A.  and  U„ 
for  all  n.  is  zero  for  S = j.  Hence,  in  (AT)  the 
cubic  and  higher-odd-order  terms  in  (4,;  are 
nonzero  only  if  s -•  1 . 

In  the  simplified  mean-field  approximation  for 


I 


41 


U BIQUADRATIC  EX C II  AN  G) 

which  the  oigencncrgicH  for  the  paramagnetic 
state  arc  given  in  Eq.  (II),  we  need  only  tie  con- 
cerned with  (S,)  as  expressed  in  (A3a),  but  where 
now 

*W4a//.W0U S,),  Y-.J'a(S,)‘.  (A8) 

Again,  lor  jJ(.Vm  +Yin~)<<  l,  Eq.  (A3a)  can  be 
written  as  (A4a),  from  which  by  virtue  of  (A8) 
we  proceed  directly  to  (A5).  The  latter  equation 
is  again  expressible  as  (A6),  from  which  we  now 
obtain 

£PflWe,i  =gHaH,+J +Pj(s«)5*  ••• . (A9) 

where  the  p’s  have  the  same  form  as  in  (A7),  but 
with  different  numerical  coellicients,  i.e., 

p,  = £a; fcwyv;,  p,  * £ b; w0)"J'a, 
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Ai«i  V,,  A'„  = 0 for  n l; 

«i  » I , A-vu/*^  * 1 >•  "1  * - 1 Av  < ''-.i/ir^'  ♦ 1 ). 

B'„  = 0 for  1*2, 

where  A,  is  defined  as  before.  Since  Av<  =0  for 
S'  -i,  we  again  have  that  in  the  expression  for 
the  cubic  and  higher-odd-order  terms  in 
($t/  are  nonvanishing  only  if  S*l. 

appendix  b 

Although  the  Hamiltonian  in  Eq.  (15)  for  the  or- 
dered magnetic  state  is  nondiagonal,  we  are  for- 
tunately interested  in  the  relatively  simple  case 
of  S=  1.  In  this  case,  the  cubic  secular  equation 
for  the  eigenenergies  is  as  follows: 


E^2Cl(SI)s*(S, 


„ - [a J<s,)2  +B\s,y  - c2((s,y 


{S.)*)1-  (C* - D3)<,S,y(S.)’]E. 


-l(A!  + B‘)C-2.\ED-IC2-D2)C({S')2  + (S.)2)\(S,)\S,y  = 0,  (Bl) 


where 

A=Jy+J2+2pBHx,{Sx),  B = J,  -J 
C=J\*J'2, 

In  principle,  with  the  solutions  for  £„,  the  magnet- 
ization can  be  obtained  from  Us  statistical-me- 
chanical definition, 

M,  = 2 pB{Sx  > = - £ (—’)  ^ (B2) 

However,  the  general  solution  for  \lx{Hx,  T)  is  a 
very  complicated  analytical  problem  and  will  not 
be  pursued  further  here. 

It  serves  our  present  purposes  to  consider  the 
situation  at  T = 0,  for  which  (B2)  reduces  to 


since  3 £,„,  3H,  = 0.  From  the  definitions  of  A and 
Mt,  the  latter  equation  can  be  written  as 


d£|  s . i c y 


from  which  we  find 


A - B *C  - D - 2(C  - D>[1  ♦ (C  - D)  2B  ](S,)2  •<■  • • • • 


Thus,  from  the  definitions  of  A,  B,  C,  and  D, 
pB  H,  - - (J2  -J'2KS,)  - 2J1(1  *■« J\ 


where  £,  is  the  ground-state  energy.  Further- 
more, at  T = 0,  we  also  have  the  simplifying  re- 
lationship 

<S„)MS,>2  = 1, 

which  we  use  to  eliminate  {Sx)  from  (Bl).  Follow- 
ing this  tactic,  we  solve  (Bl)  and  obtain  £,  as  a 
power  series  in  (Sr)2,  namely, 

£1  = £,0  + £lJ(S,)a^£1<(SI)V...,  (B4) 

where  £I0  = -B-C  and  Eu,  £,„•••  are  compli- 
cated functions  of  A,  B,  C,  and  D Substitution 
of  (B4)  into  (B3)  gives 


which  leads  directly  to  Eq.  (16)  for  Mt  vs  Hx  at  T 
= 0. 

In  the  absence  of  Hx,  (5,)  = 0 and  (Bl)  reduces  to 

El  * 2 c(s,)2e2„  - U?\s,>2  - c2(  s,y]  £„  = 0, 

whose  solutions  at  e 

£„  = 0,  ±B\SX)  - C\S,)2.  (B5) 

Furthermore,  since  the  Hamiltonian  (15)  is  now 
diagonal  in  ( Sx),  it  follows  from  Eqs.  (B5)  and 
(6a)  that 

#,BB(Sf>  _ |,-«8(  Sg) 

(S,)*eaB(S.>  +e-IC(St>J  (D6) 

which  gives  the  sublattice  magnetization  (2 pa\S,)) 


as  a function  Of  temperature  up  to  the  N£el  point 
(Ty).  To  obtain  T v,  we  let  ( S, ) in  (D6)  beconie 
vanishingly  small  and  find  ihat 

T s = 2B/2k  = 2(«/,  2k.  (07) 

For  the  susceptibility  at  T.v,  we  let  (S,)  = 0 m 
(Bl),  which  now  reduces  to 

El*2CiSsyEi„-[A^St)1-CiiS,y)E„  = 0. 

Hence,  the  eigenenergies  are 

£,  = 0,  ±A(Sz)-C{S,y, 
which,  substituted  into  (B2),  yield 


ctf,MSx>  _ 

V ♦<r4v‘,lV  ,t,-avc<st>-  ’ (DtJ) 

where  3V  = 1 kT  s.  Through  the  dependence  of  A on 
//,,  (B8)  describes  the  induced  magnetization 
(2 M8(S,))  at  7,  as  a function  of  //,  For  very 
small  values  of  //,  and  \S, ),  (B8)  reduces  io 

(Sl)  = 5^A<SI)*2[(./,WIKS,)  <•  2/ifl //,  j 3 AT., 

which,  combined  with  (B7),  gives 

X~2(iat,Sx),  H r = — 2pa  V2 

for  the  initial  susceptibility  at  T s. 
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Appendix  F 

FIELD-INDUCED  TRANSITIONS  IN  DySb* 


T.O.  Brun,  C.H.  Lander,  F.W.  Korty  and  J.S.  Kouvel 
Argo one  National  Laboratory,  Argoone,  Illinois  60439  and 
University  of  Illinois,  Chicago,  Illinois  60630 


ABSTRACT 

The  NaCl-structured  compound  DySb,  which  In.  zero 
field  transforms  abruptly  at  Tjj  % 9.5  E to  a Type- II 
anciferromagmetic  (A)  state  with  a nearly  tetragonal 
lattice  distortion,  was  previously  found  to  exhibit 
rapid  field- Induced  changes  in  magnetization  at  1.5  E. 
The  field-induced  transitions  in  a DySb  crystal  have 
been  studied  by  neutron  diffraction  and  magnetization 
measurements  in  fields  up'  to  i 60  kOe  applied  parallel 
to  each  of  the  principal  axes.  In  the  <100>  case,  the 
transition  from  the  A to  an  intermediate  ferrlsagnetlc 
(Q)  state  is  first-order  at  4.2  K (critical  field  Hc  % 
21  kOe)  but  is  continuous  from  h,  6 K up  to  T^,  as 
Ke  — 0.  The  Q-to-paramagnstlc  (?)  transition  1s  rapid 
but  continuous  at  4.2  & (Hc  i 40  kOe)  and  becomes 
broad  as  T«j  is  approached.  In  the  <110>  case  the  A-Co- 
Q transition  remains  essentially  first-order  from  4.2 
K (Hc  q,  15  kOe)  up  to  Ts;  above  Tjj  rapid  P-co-Q  tran- 
sitions occur  at  very  high  fields.  The  magnetic  struc- 
ture of  the  Q state  is  found  to  be  that  of  SoP. 

INTRODUCTION 


In  zero  magnetic  field,  the  compound  DySb 
exhibits  a first  order  megnetlc  transition-'^  at  T^  % 
9.5  K.  Above  T«j  DySb  has  the  cubic  NaCl-structure. 
Below  Che  compound  orders  anciferromagnecically  in 
the  type-II  structure  that  consists  of  ferromagnetic 
(111)  planea  stacked  antif erromagaeClcaUy  along  the 
[111]  axis.  The  first  order  transition  is  accompanied 
by  a predominantly  tetragonal  lattice  distortion  (e/a  ■ 
.993).  The  ordered  magneclc  moment  at  6 E la  9.5  Un 
per  Dy,  almost  the  saturation  momant  of  10  Ug  for  Dy  , 
and  Is  parallel  to  the  tetragonal  [001]  axis.  Elastic 
as  well  as  aagneto-tharmal  data  for  DySb  have  been 
analyzed  in  a molecular  field  approximation, and 
tha  analysis  shows  that  biquadratic  pair  interactions 
between  tha  Dy-1+  loos  may  ba  as  important  as  the 
bilinear  exchange. 

Previous  magnetization  oaasurements^  made  on  a 
single  crystal  of  DySb  ac  1.5  E in  fields  up  to  60 
SSi  show  that  * ferriaagnetic  state  exists.  For  all 
thrse  principal  directions  of  tha  field  Che  values  of 
Che  aagnacizatlon  in  this  state  are  consistent 
with  Che  megnedc  structure  of  DySb  being  Chat  of 
BoP,^  in  which  the  moments  within  each  (111)  plana  of 
Dy  atoms  are  ferromagnatlcally  aligned  with  10  u3  per 
Dy  and  osclllaca  between  two  perpendicular  <100> 
directions  In  alternating  (111)  p lanes . 


RESULTS 


Magnetization  measurements  were  performed  on 
elaglc  cryscals  of  DySb  using  a vlbratlng-sample 
magnetometer.'  Magnetic  fields  of  up  to  56  kOe  were 
applied  parallel  to  tha  <100>,  <110>  and  <111> 
crystallographic  directions.  The  fields  have  been 
corrected  for  demagnetization.  The  results  for  tha 
<100>  direction  are  shown  in  Fig.  1 as  iaocharms  M vs 
H.  Below  T^  we  observe  two  critical  flalds  at  which 
the  sagaatlzaclon  increases  very  rapidly.  Ac  4,2  £ the 
magnetization  changes  discontinuous  ly  at  Hc,<1,00>  - 
21.6  kfle  and  very  rapidly  at  “ 40.6  kOe. 

Both  critical  fields  dacresse  with  increasing  tempera- 
ture, both  transitions  becoming  broader  as  the 
temperature  approaches  T^.  At  4.2  K the  magnetization 
is  almost  5 Ug  per  Dy  atom  Just  above  R?i<l“0>  end 
approaches  10  Ug  par  Dy  atom  above  Hgj**™*. 

In  Fig.  2 the  M vs  B isotherms  are  shown  for  the 
field  along  <U0>.  The  iaocharms  show  one  transition 
which  remains  abrupt  while  lea  critical  field, 
aci<110>,  decreases  from  15.6  kOe  at  4.2  K to  zero  at 
T,;.  The  magneclzaclon  at  4.2  K approaches  e value  of 
approximately  7.2  Ug  at  high  fields.  In  addition,  the 


H (kOe) 

Figure  1.  Magnetization  of  DySb  vs.  field  applied 
along  <100>  ac  constant  cemperacure. 


Figure  2.  Same  ea  Fig.  1 for  field  along  <100>. 

isotherm  for  9 K shows  e very  rapid  Increase  In  CJ  d2 
at  flelda  Juat  above  • The  10  E and  12  E iso- 

therms above  T^  show  a similar  change  In  dM/dH  ac  flelda 
Just  above  15"  kOe  and  26  kOe,  respectively. 

To  determine  the  ranges  of  temperature  and  field 
for  the  various  magneclc  states  of  DySb,  we  have  con- 
verted our  data  to  field  vs  temperature  curves  of 
constant  aagnacizatlon.  The  phase  boundaries  have  beta 
Identified  from  the  rapid  changes  in  dH/dT.  Fig.  3 
shows  the  magnetic  phase  diagram  for  H parallel  to 
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2<t$ 


Figure  3.  Field  along  <100>  vs.  temperature  at 
constant  magnetization  (in  Pfl/Dy)  for  DvSb. 


Figure  A.  Same  as  Fig.  3 for  field  along  <110>. 
<100>.  The  phase  marked  P Is  paramagnetic,  the  A phase 
is  the  Type- II  antiferromagnetic  state,  and  the  Q-phass 
is  the  intermediate  ferrlaagnetlc  stjte. 

The  magnetic  phase  diagram  for  H parallel  to  <110> 
is  shown  in  Jig.  A.  The  P,  A,  and  Q phases  are  the 
same  ea  for  B parallel  to  <100>.  An  unusual  feature  of 
this  diagram  is  that  the  boundary  between  the  Q and  P 
phases  extends  up  to  15  K at  the  maximum  field. 

The  magnetization  data  for  H||<111>  are  not  shown, 
but  the  magnetic  phase  diagram  is  similar  to  that  for 


H||<110>.  The  critical  field,  Hcl<m>,  is  19  kOe  at 
A. 2 K,  and  the  value  of  the  magnetization  in  the  Q 
phase  is  n,  5.8  per  Dy  atom. 

Neutron  diffraction  measurements  have  been  per- 
formed on  a single  crystal  of  DySb  with  ff||  [110]  to 
determine  the  structure  of  the  Q-phase.  We  have 
measured  the  integrated  intensities  of  the  (y  7 7) 
and  the  (5^1)  antiferromagnetic  reflections  at  4.2  K 
and  11  K ana  for  H “ 0 and  30  kOe.  The  observed  in- 
tensities for  the  Q-phaae  are  in  good  agreement  with 
10  Pn  per  Dy  parallel  to  [100]  and  [010]  in  alternate 
(111)  planes.  At  11  K no  magnetic  reflections  are 
present  In  zero  field,  but  by  increasing  the  field  the 
(i  7 7)— reflection  appears  at  the  P to  Q boundary, 
with  the  same  Intensity  as  at  4.2  K.  Preliminary  data 
for  the  lattice  distortions  in  the  Q phase  indicate 
that  the  lattice  expands  in  the  [001]  direction, 
normal  to  the  plane  of  the  moments. 

DISCUSSION 


The  results  of  the  present  magnetization  measure- 
ments are  in  good  agreement  with  the  results  of  Busch 
and  Vogt.^  As  pointed  out  by  these  authors  the  value 
of  the  magnetization  in  the  Q phase  for  each  of  the 
three  principal  field  directions  is  consistent  with 
the  HoP  structure.  Furthermore  if  the  three  Q phases 
have  identical  magnetic  structures  with  the  same  ex- 
change and  anisotropy  energies,  the  differences  be- 
tween the  Zeeman  energies  of  the  three  Q-phases  pre- 
dict that  H fl<100>  “ ^ Hcl<110>  - 2/d'H^1^.  The 
experimental  values  of  the  three  critical  fields  are  in 
good  agreement  with  this  prediction.  Both  magnetiza- 
tion and  neutron  diffraction  results  for  the  Q-phase 
of  DySb  show  therefore  that  this  phase  has  the  HoP 
structure  with  10  per  Dy  atom. 

For  J||<100>  the  HoP  structure  cannot  be 
stabilized  by  bilinear  exchange  Interactions  alone; 
additional  couplings  such  as  biquadratic  pair  inter- 
actions are  needed.  Stevens  and  Pytte^  have  suggested 
that  the  dominant  contribution  to  this  biquadratic 
interaction  comes  from  the  couplings  between  the  Ion 
and  the  lattice.  In  cheir  model  cooperative  shifts 
of  the  Sb  ions  relative  to  the  Dy  ions  create  local 
distortions  stabilizing  the  HoP  structure. 

Further  neutron  diffraction  measurements  are  re- 
quired to  verify  the  HoP  structure  of  the  Q phases  for 
other  direction  of  the  field  and  to  test  the  model  of 
Stevens  and  Pytte?  Our  magnetization  data  for  DySb 
above  TN  are  being  compared  in  detail  with  crystal- 
field  calculations  of  magnetization  vs  effective 
field  in  order  to  expose  any  evidence  for  biquadratic 
or  any  other  higher-order  coupling,  as  was  done  In  a 
recent  investigation  of  PrAg. 

*Based  on  work  performed  under  the  auspices  of  the 
U.  S.  Atomic  Energy  Commission. 

‘''Supported  In  part  by  the  Office  of  Naval  Research. 
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BIQUADRATIC  EXCHANGE  INTERACT!'  INS  IN  DySb* 

J.S.  KOUVEL.+  T.O.  BRUN  and  F.W.  KORTY  + 

Argonne  Salional  Laboratory,  Argonne.  Illinois  AfW.tV.  USA  and  Unit-mils  of  Illinois.  Clin  ago,  Illinois  hitbSn.  USA 


From  a comparison  of  magnetisation  data  on  a DySb  crystal  above  its  Neel  point  with  magnetisation  calculations 
based  on  its  crystal-field  states,  it  is  shown  that  the  exchange  interactions  have  a substantial  biquadratic  component,  as 
suggested  by  the  field-induced  stale  with  orthogonal  sublattice  moments. 


At  the  Neel  temperature  (TS=9.5K)  in  zero 
magnetic  field,  the  cubic  NaCl-structured  com- 
pound DySb  transforms  abruptly  into  a type-ll 
antiferromagnet  in  which  the  Dy  moments  in 
adjacent  ferromagnetic  (III)  planes  are  an- 
tiparallel (I,  2].  The  first-order  transition  is  ac- 
companied by  a nearly  tetragonal  lattice  distor- 
tion and  the  ordered  moments  (~10aib/Dv- 
atom)  lie  along  the  tetragonal  [0011  axis. 
Moreover,  magnetic  measurements  [?]  on  a 
DySb  crystal  at  1.5  K have  shown  that  the 
compound  undergoes  a rapid  field-induced 
transition  from  the  antiferromagnetic  state  to  a 
state  whose  magnetization  corresponds  to  an 
alignment  of  Dy  moments  in  equal  numbers 
along  each  of  two  orthogonal  (100)  directions.  It 
was  suggested  [ 3 1 that  the  magnetic  structure  of 
this  quadrature-spin  (Q)  state  of  DySb  was  that 
of  HoP,  in  which  the  moments  in  adjacent  fer- 
romagnetic (111)  planes  are  aligned  ortho- 
gonally, and  this  has  been  borne  out  by  neutron 
diffraction  measurements  on  a DySb  crystal  in 
high  fields  14 J.  The  latter  work  |4|  also  included 
a magnetization  study  which  gave  a deter- 
mination of  the  magnetic  phase  diagrams  for  the 
applied  field  I H)  along  < 1 (K)>  and  <1 10).  These  are 
shown  in  fig.  1.  together  with  the  phase  diagram 
for  H along  (111),  which  we  recently  deter- 
mined from  a similar  study  of  the  same  DySb 
crystal.  The  nature  and  positions  of  the  phase 
boundaries,  plus  the  occurrence  of  tricritical 
points  (51,  will  be  discussed  in  a separate  report. 

For  our  present  purpose,  the  magnetic  phase 
diagrams  for  DySb  serve  to  emphasize  the  high- 
field  existence  of  the  Q state,  which  cannot  be 
stabilized  by  bilinear  exchange  interactions 
alone  but  requires  higher-order  interactions  in 

•Work  performed  under  the  auspices  of  ihe  US  Energy 
Research  and  Development  Administration 

♦Supported  in  part  by  the  National  Science  Foundation 
and  the  Office  of  Naval  Research. 


<100>  <110>  <111> 


T (°K) 

Fig.  1 Magnetic  phase  diagrams  of  field  vs  temperature  for 
DySb.  indicating  the  anliferromagnenc  lA).  quadrature-spin 
(Ql.  and  paramagnetic  (P)  regimes. 

addition,  as  was  pointed  out  earlier  |4|.  Spe- 
cifically, if  the  higher-order  interactions  have  a 
simple  biquadratic  form,  the  exchange  hamil- 
tonian  (of  the  ith  spin)  can  be  written  as 

K = - 2 I-7-5.  • s,  +• j is,  ■ s,  ri.  ( 1 ) 

I 

The  existence  of  the  Q state  would  then  suggest 
that  the  inter-sublattice  component  of  the 
biquadratic  exchange  coefficient  J'„  is  negative, 
thus  favoring  quadrature  (orthogonal)  over  col- 
Itnear  alignment  of  the  sublattice  moments.  The 
importance  of  biquadratic  interactions  in  DySb 
has  also  been  deduced  previously  from  a mean- 
field  analysis  of  elastic  and  magnetothermal 
data,  and  it  was  concluded  that  J\,  is  effectively 
positive  (6,  7).  Although  a positive  provides  a 
plausible  rationale  for  the  first-order  transition 
of  DySb  at  TN  in  zero  field  (7),  it  was  shown 
more  recently  that  the  abruptness  of  this  tran- 
sition can  be  justified  solely  on  the  basis  of 
magnetic  symmetry  [8). 
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Fig  2.  Calculated  (dashed)  curses  of  Vf!  vs  H.,IM  and 
experimental  data  of  ,\f!  vs  HIM  for  DvSh  at  IS  K svilh  the 
effective  field  (H„ I.  applied  field  tH).  and  magnetisation 
l At  I along  various  crystallographic  directions 


In  order  to  reveal  explicitly  any  biquadratic 
coupling  in  DvSb.  we  have  followed  a pro- 
cedure recently  used  for  the  compound  PrAg 
|9],  First,  we  had  to  determine  the  mag- 
netization ( M ) vs.  effective  field  tH,„)  behavior 
of  DySb  from  its  crystal-field  slates.  The  latter 
are  presently  available  from  interpolations 
between  the  neutron  scattering  results  on 
various  other  isomorphic  rare-earth  antimonides 
( 10).  We  could  thus  establish  quite  reliably  that 
the  crystal-field  coefficients  for  DySb  are 
A"(r‘)  = 60  K and  A»(r6)  = 2 K,  from  which  we 
deduced  that  relative  to  the  lowest  lying 
doublet  the  three  quartets  are  at  8.8.  96.0  and 
120.2  K.  and  the  F-  doublet  is  at  60.4  K [ll|. 
From  this  crystal-field  scheme,  we  calculated  .V/ 
vs.  H,«  at  various  temperatures  for  //<.„  along 
each  of  the  principal  crystallographic  directions, 
the  results  for  18  K are  presented  as  M'  vs. 
H'tIM  in  fig.  2.  Also  shown  in  this  figure  are  our 
experimental  magnetization  results  for  DySb  at 
18  K in  fields  up  to  56  kOe  applied  along  the 
same  three  directions,  which  are  plotted  analo- 
gously as  iVf:  vs.  HIM.  These  and  similar  ex- 
perimental curves  obtained  at  other  tem- 
peratures within  the  paramagnetic  regime  were 
each  compared  against  the  corresponding  cal- 
culated ,\f;  vs.  H'tlM  isotherm; 
was  then  evaluated  at  each  experimental  value 
of  M2. 

Our  results  are  plotted  as  M:  vs.  ( H,*  - H)IM 
in  fig.  3.  They  can  be  compared  meaningfully 


(HetrH)/M  IkCe/ag'Dy) 

Fig  l \f  vs.  (H,«-H)/.\f  for  DvSh  at  various  tem- 
peratures (in  Ki  for  different  directions  of  the  effective  field 
i H.«».  applied  field  (H  I.  and  magnetization  <M  i 

with  an  expression  for  the  total  exchange  field, 

H„ctl  = - H = A\f  + A ' Af  ’.  (2) 

which  has  been  derived  |9|  from  spin  Hamil- 
tonian (I)  within  a mean-field  framework, 
where  A and  A’  are  proportional  to  the  /-sum- 
mations of  J„  and  J\.  respectively,  for  the 
paramagnetic  state.  Thus,  if  A and  A'  were 
constant.  ,Vf:  vs.  (Hrlf-H)/M  would  describe  a 
straight  line  and,  if  A'  were  zero,  the  line  would 
be  vertical.  The  pronounced  deviations  of  the 
curves  in  fig.  3 from  verticality  can  therefore  be 
taken  as  evidence  for  substantial  biquadratic 
interactions  in  DySb.  Furthermore,  these  in- 
teractions vary  significantly  with  temperature 
and  with  the  direction  (and.  in  some  cases,  the 
magnitude)  of  the  magnetization.  A detailed 
analysis  of  these  results  will  be  presented  in  a 
future  report. 
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